LU TP 06-16 
|hep-ph/0604043| 

revised July 2006 



O : Chiral Perturbation Theory Beyond One Loop 

(N 

wo: 
< 

Johan Bijnens 

(N 

>. 
CO . Department of Theoretical Physics, Lund University 

^ ! Solvegatan 14A, SE 22362 Lund, Sweden 

O 

\o 
o 

22 

t^ : 

> 

X 

j^ ■ Abstract 

The existing Chiral Perturbation Theory (ChPT) calculations at order p e are reviewed. The 
principles of ChPT and how they are used are introduced. The main part is a review of the two- 
and three-flavour full two-loop calculations and their comparison with experiment. We restrict 
the discussion to the mesonic purely strong and semileptonic sector. The review concludes by 
mentioning the existing results in finite volume, finite temperature and partially quenched ChPT. 



Chiral Perturbation Theory Beyond One Loop 

Johan Bijnens 

Department of Theoretical Physics, Lund University 
Solvegatan 14A, SE 22362 Lund, Sweden 

February 7, 2008 

Abstract 

The existing Chiral Perturbation Theory (ChPT) calculations at order p e are reviewed. The 
principles of ChPT and how they are used are introduced. The main part is a review of the two- 
and three- flavour full two-loop calculations and their comparison with experiment. We restrict 
the discussion to the mesonic purely strong and semileptonic sector. The review concludes by 
mentioning the existing results in finite volume, finite temperature and partially quenched ChPT. 

1 Introduction 

Chiral Perturbation Theory (ChPT) is the low-energy effective field theory of Quantum Chromo Dy- 
namics (QCD) where the degrees of freedom taken into account are the Goldstone bosons from the 
spontaneous breakdown of the chiral symmetry and their interactions. The subject grew out of the 
current algebra approach of the 1960ies. It was brought into its modern form by Weinberg, Gasser and 
Leutwyler PUEHHJ- Especially the work of Gasser and Leutwyler was instrumental in the renaissance 
of effective field theory methods in low-energy hadronic physics. Introductions to ChPT can be found 
in the lectures of Refs. jHEHH] as weu as m Sect. 121 Shorter introductions can be found in Refs. [ZUHJ. 
The lectures by Leutwyler [5| stress the foundational aspects of ChPT and are a very recommended 
read. Introductions to QCD with an emphasis on the low-energy aspects and ChPT can also be found 
in the books by Donoghue, Golowich and Holstein [TO], Georgi [TT] and Smilga [T2] . 

Chiral Perturbation Theory (ChPT) is now a very large subject, no single review can do the entire 
field justice. Two main earlier reviews when the main part of ChPT was completed to next-to-leading- 
order (NLO) or order p 4 are those by Ecker J3] and MeiBner |14j . This review concentrates on results 
at next-to-next-to-leading-order (NNLO) or order p 6 in mesonic ChPT, restricted to purely strong 
or electromagnetic processes. The wide field of applications beyond this, weak mesonic decays JEHU], 
electromagnetic corrections to hadronic processes [To^ lH)]. non-relativistic approaches to hadronic atoms 
[TTj and the entire area of baryon and nuclear physics ^H] are not treated. The articles cited present 
an entrance to those fields. This review includes the comparison with and prediction of experimental 
data in connection with the existing order p 6 calculations. 

In Sect. |2] ChPT itself is discussed. I start there by discussing the global symmetries of QCD 
in Sect. 12.11 The consequences of a global symmetry are embodied in the Ward identities involving 
the Green functions of the theory. A very elegant method to automatically produce Green functions 
obeying the Ward identities is the external field method discussed in Sect. 12.21 In QCD, the chiral 



symmetry is present in the Lagrangian but it is not visible in the spectrum. The way to combine those 
two observations is by looking at the concept of spontaneous symmetry breaking introduced here in 
Sect. 12.31 In fact, the chiral symmetry is not quite exact in QCD. It is only valid when all quark masses 
are zero. But for the light quarks, this can be treated as a perturbation. The way of dealing with this 
in general is discussed in Sect. 12 .41 and the spontaneous symmetry breaking in QCD is briefly discussed 
in Sect. 12.51 After this I proceed with ChPT at lowest order, built up from the Goldstone bosons from 
spontaneous chiral symmetry breakdown in QCD. The lowest order is introduced in Sect. 12.61 and its 
main applications at that order are shown in Sect. 12.71 Nonrenormalizable effective field theories, as 
ChPT, have in principle an infinite number of parameters. In order for them to be phenomenologically 
useful, there has to be a way to order the various parts in order of importance. This concept is called 
powercounting and discussed in Sect. 12.81 Its interplay with renormalization is also introduced there. 
Once the concept of powercounting is established, one has to find out how to construct the explicit 
Lagrangians needed at each order in the powercounting and find a way to do the renormalization in 
practice. The main ideas behind both those subjects are treated in Sects. 12 . 91 and 12 . 1 (Jl respectively 

After this brief introduction to ChPT and the principles behind it, the main part of this review 
follows. It is split into three parts. First a review of the existing calculations in infinite volume for the 
case of two flavours. Here only the up and down quark are treated as light and the relevant degrees 
of freedom are the pions only. The next part treats the three-flavour case. Here up, down and strange 
quarks are all treated as light and the full lowest mass pseudoscalar octet of pions, kaons and eta are 
included as the relevant degrees of freedom. This part includes an overview of all order p 6 calculations 
at infinite volume relevant in this domain. The third part consists of those order p e calculations which 
do not belong in either of the two previous categories. It includes finite volume and finite temperature 
calculations. I also briefly discuss what is known for the partially quenched regime there. 

For the two-flavour case I first discuss the early estimates of NNLO effects using dispersive methods 
in Sect. 13.11 The remaining subsections go through the various processes known to NNLO order. 
77 — ► 7T°7r was the first case of a full NNLO calculation and is discussed in Sect. 13.21 The pion mass 
and decay constant follow in Sect. 13.31 77 — » 7r + 7r~ and the pion polarizabilities are an area of active 
experimental work and at present there seems to be a discrepancy between data and ChPT. The relevant 
calculations are reviewed in Sect. 13.41 A major recent success of theory of low-energy hadronic physics 
is the extremely accurate description of low-energy pion-pion scattering. The role of ChPT and the 
relevant calculations are discussed in Sect. 13.51 We have also included the full NNLO ChPT formula 
here since it can be beautifully expressed in terms of elementary functions. The two remaining existing 
calculations are those of the pion vector and scalar form- factors, reviewed in Sect. 13.61 and the pion 
radiative beta decay, n — >■ £uj. The latter is discussed in Sect. 13.71 The present best values of the 
low-energy constants (LECs) at NNLO are given in Sect. 13.81 

When one considers three light quark flavours, the number of processes increases rather dramatically. 
Due to the presence of many different scales, the calculational difficulty also increases. Nonetheless, a 
great many processes are known to NNLO also in this sector. This starts with the vector two-point 
Green functions, Sect. 14.11 which were used as a laboratory for doing three-flavour NNLO calculations 
and some small phenomenological applications. A very similar calculation is required for the scalar 
two-point functions discussed in Sect. 14.21 The simplest two- loop calculation is the quark- ant iquark 
condensate. It is treated in Sect. 14.31 The first calculations in the three-flavour sector requiring proper 
or irreducible two-loop integrals are the axial- vector two-point functions, the pseudoscalar meson masses 
and the decay constants. These are reviewed in Sect. 14.41 The next calculation is in fact one of the 
more elaborate ones, for the decays K — > kkIv. The reason this was done was that this process is one 
of the main sources of the order p A or NLO LECs. Results are reviewed in Sect. 14.51 With those results 
in hand, the basic set of parameters of ChPT could be determined to NNLO. This lead to a full flurry 
of applications. All electromagnetic, Sect. 14.61 and scalar form-factors, Sect. 14.81 of the pseudoscalar 
mesons have been worked out and compared with existing data. The process K — > nlv contains several 



form-factors and is needed for the determination of the CKM-matrix element V us . The ChPT aspects 
of this are reviewed in Sect. 14.71 This section proceeds with what is known for pion-pion scattering and 
pion-kaon scattering and some possible consequences for the LECs L\ and L r % . As very briefly indicated 
there, this is relevant for a possible strong flavour dependence of spontaneous chiral symmetry breaking. 
We present some results here as well from the last remaining calculation, ir, K — > ivy, in Sect. 14.111 
We conclude with a few comments about the estimates of the NNLO or order p 6 LECs. This is one of 
the main open questions in this field. Note that there have been a few claims of two-loop calculations 
in the literature which neglected the proper two-loop diagrams. None of these are mentioned in this 
review. 

The remaining part of this review is kept very short. I basically only mention the calculations which 
have been done for finite temperature and volume and give only a brief overview of the order p 6 work 
done for the partially quenched case. 

A few last remarks, a website with links to the actual formulas of many of the papers reviewed here 
as well as some lectures on ChPT and more general effective field theory is Ref. [T^]. There are also 
many calculations in the anomalous sector. Here the order p 6 is only one-loop. A review is Ref. [2U| 
and references to more recent work can be found in the papers place where the order p 6 Lagrangian 
for this sector was worked out [2U 122] • I have done a reasonable effort to dig out all relevant papers 
for this review. With ChPT being such a large subject, I have however most likely overlooked some 
directly relevant for the subject considered. 

2 Chiral Perturbation Theory 

2.1 Chiral Symmetry 

At the low energies discussed in this review only the lightest quarks are relevant. The heavier quarks, 
charm, bottom and top play no role here. We put the quarks together in a column vector 




9=\ , ! (1! 

for the two-flavour case or 

for the three-flavour case. The conjugate row vector is analogously defined via 

q=(ud) or q=(uds). (3) 

The generalization to n F flavours is obvious. 

The gluons couple identically to all quark flavours. If all the masses are equal for n F flavours we 
have a SU{n F ) v x U(l)y symmetry. The U(l)y symmetry changes the phase of all the quark fields 
simultaneously and corresponds to baryon number. The SU(n F )y symmetry acts as 

q(x) ► Uq(x) U G SU(n F ) v . (4) 

The vector symmetry SU(n F )y is known as isospin for the two-flavour case and as the Gell-Mann- 
Ne'eman octet symmetry for the three-flavour case. 

However, QCD has a larger symmetry structure. The QCD Lagrangian is of the form 

£qcd = J2 tQiLpQiL + iQiR^QiR ~ miq iR q iL - miq iL q iR + ■■■ . (5) 

i=u,d,s 



Here Jp is the covariant derivative with the gluon field and the dots indicate the purely gluonic terms. 
The sums over colours are understood and not explicitly written out. The left and right handed quark 
fields are given by 

Qr = ~ (1 + 75) q and q L = - (1 - 75) q . (6) 

We define a quark mass matrix 

j m u \ 

M=\ m d \ (7) 

V m s J 

and left and right-handed column vectors q^ and q R by replacing in (JTJ) and (J2J) all quark fields by their 
left and right-handed parts. This allows us to rewrite the QCD Lagrangian as 

£qcd = HiJP Ql + iq R Ip qn - q R Mq L - q L Mq R . (8) 

The form (jHJ) shows that there is in fact a larger symmetry than the flavour rotations of (JIJ) whenever 
the quark masses are equal to zero, M. = 0. To be precise, we obtain the chiral symmetry group 

G = SU(3) L x SU(3) R x U(l) v x U(l) A , (9) 

with 

qR — > 9RqR 9r e SU(3) R 

q L -^g L q L g L eSU(3) L . (10) 

Under U(l)v all quarks have the same change in phase while under U(1)a the right and left-handed 
quarks have the opposite change in phase. The symmetry is called chiral because it acts differently on 
the left and right-handed quarks. 

The U(1)a is only a symmetry of the classical action, not of the full quantum theory of QCD. The 
divergence of the associated current does not vanish due to the anomaly |2SJ- It is nonzero by a total 
divergence but instantons allow for this to have a physical effect [23J- We will not consider its effects 
for the remainder of this paper. 

The U(l)y symmetry corresponds to baryon number. We will only discuss mesons in this review 
and hence we also drop this symmetry. The final chiral symmetry of QCD in the limit where all quarks 
are massless is thus 

G x = SU(3) L xSU(3) R . (11) 

2.2 External field method 

The consequences of a global symmetry are most clearly expressed through the Ward identities for Green 
functions. These can be derived using the methods given in most field theory books but a particularly 
elegant method to obtain Green functions that obey the Ward identities is the external field method. 
This method also allows to show clearly how the knowledge of the Green functions of QCD in the chiral 
limit is sufficient also to describe the Green functions away from the chiral limit. The particular version 
described here was introduced by Gasser and Leutwyler [2]. 

For the mesonic physics we discuss in this paper we look at Green functions or correlation functions 
defined by vector and axial-vector currents, and scalar and pseudoscalar densities. These are referred 
to together as external currents. The currents are defined by 

v?(x) = q^h^qjix) , 
A% ii(. x ) = Qi(xhtf5qj(x) , 

S l3 (x) = -^(x)^-(ar) , 

P«(x) = q^ij^ix) . (12) 



The indices i,j run over the quark flavours u, d or u, d, s and the sum over colours is implicitly under- 
stood. 

Green functions or correlation functions can in general be introduced via the inclusion of sources 
in the Lagrangian. This is described in most books on Quantum Field Theory, see e.g. chapter 9 in 
Ref. |25j . These sources are referred to as external sources or external fields. 

The Lagrangian of massless QCD extended by the external fields is written as 

£qcd = nJP Ql + iq R Ip qR-q R (s + ip)q L -q L (s- ip) q R + q L ^l^q L + q R l^r^q R . (13) 

The external fields s, p, l^ and r^ are space-time dependent np x rip matrix functions. The vector and 
axial- vector fields, v^ and a M , are included via 

^ = f M -o M , r^^v^ + a^. (14) 

These external fields are all Hermitian matrices. 

By taking functional derivatives with respect to the external sources s, p, v^ and a^ from the 
functional integral we can then build up all the wanted Green functions with the currents defined 
in ()12|) of massless QCD. To get an insertion of the current V^(x) one needs to take the functional 
derivative w.r.t. v^{x). 

With these additional sources the global chiral symmetry described in the previous section can be 
extended into a local chiral symmetry. A local symmetry transformation is given by an element of the 
symmetry group, gp x g R G SU(tif)l x SU(np) R where gi and g R are now functions of the space-time 
point x. The transformations under the local symmetry are 



Ql - 


-> 9hqi 


qn - 


■> 9B.qR 


\4 = (s + ip) - 


■* 9RMg[, 


'/i = V/J, Cifl ~ 


-> 9Ll^9L- id ^9L9L^ 


r^ = v^ + a^ - 


"> 9Rr^g R -id^g R g R . 



(15) 

For most of this paper, we will consider these symmetries as exact also at the quantum level. They are 
anomalous but that effect is fully taken into account by the Wess-Zumino-Witten|2Ell2Zj term and that 
one can be taken explicitly into account as described in [2J E] . 

The function that gives all the Green functions directly when the functional derivative w.r.t. to the 
external fields is taken, is called the generating functional, G and it is given by the functional integral 

\dqdqdG]e J " 
^,r»,s, P) = . ... . -<=- - x • (16) 

[dqdqdG]e ' % QCD 

Z M ,r>,s,p=0 

The function r(/ M , r M , s,p) is called the effective action. J[dqdqdG] indicates the functional or Feynman 
path integral over all possible quark, anti-quark and gluon paths or configurations. 

With this one can see how the generating functional for QCD with nonzero masses is related to the 
one in the chiral limit. We have indicated the two different cases here by the superscript m q ^ and 
m q = 0. Comparing the Lagrangians with and without the quark masses leads to the relation 
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In a similar fashion the response to external electroweak vector fields can be included 1 . The couplings 
of photons to quarks is described by a term in the Lagrangian of the type 

eA„ (q L Q L ^q L + q R Q R Yq R ) , (18) 

with 

/2/3 \ 

Ql = Qr = Q=\ -1/3 . (19) 

V -1/3 J 

A^ is the photon field. The interactions with photons can thus be included by changing 

l„ -> l„ + eA^Q , r M ^ 7> + eA^Q . (20) 

This should be understood in a way similar to the inclusion of quark masses by changing s to s + M. 
as done via Eq. (J17JI . In particular, the Green functions in the presence of external electromagnetism 
and quark masses are up to a normalization given by G(/ M + eA^Q, r M + eA^Q, s + A4,p). 

The couplings to Wr and Z^ bosons can be included in a similar fashion by looking at the standard 
model Lagrangian and seeing how those couplings can formally be written as parts of l^ and r M , just as 
we could write M. formally as a part of s. To be precise, in terms of the Cabibbo angle, 9 C , the charged 
W* 1 couplings are include by changing 



cos 6 C W+ sin0 c W+ 

J„ - J„ - -i^ | coaOcW- | (21) 

sin 6 C W^ 



9_ 
V2 



with g = e/ sin 9\y, the Weinberg angle. 

As described in this subsection, currents, external sources, external fields, electroweak gauge bosons 
are different objects but are treated in the end via the external fields l^ and r M . Many papers tend to 
be somewhat cavalier with the choice of words. 

We have described the Green functions here as calculated by using functional derivatives. This 
is of course completely equivalent to calculating them directly using Feynman diagrams. Similarly, 
amplitudes calculated directly with Feynman diagrams are equivalent to those calculated from the 
Green functions using standard LSZ reduction. The main reason for using this external field formalism 
is that it allows for calculations maintaining chiral symmetry throughout the entire calculation and 
allows for a simpler way of classifying terms in ChPT as described below. 

2.3 Spontaneous symmetry breaking 

In Sect. 12 .11 we introduced the chiral symmetry of QCD in the massless limit. This symmetry is however 
clearly not present in the spectrum of hadronic states. If it was a good symmetry realized in the same 
way as Lorentz or rotation symmetry the spectrum of strongly interacting particles would look quite 
different. In particular, one would have parity doublets. For every particle with a given spin and parity 
there would be another one with the same spin but the opposite parity. In the presence of small quark 
masses we would still expect this to be approximately true, just like we find isospin doublets and triplets 
as well SU{2>) octets and decuplets. If one looks at the mass spectrum this is clearly not the case. There 
is obviously no partner with a mass close to the pion, the rho or the proton. The possible candidates, 
ao(980), Oi and the Roper are very far away in mass and have in general very different properties. 
Clearly, the chiral symmetry, G x , is not directly visible in nature. 

^^We use here the word external since when the exchanges of electroweak bosons between strongly interacting particles 
is needed, the formalism needs to be extended. 





Figure 1: The potential V((f>) for an unbroken Figure 2: The potential V(<p) for a sponta- 
symmetry. neously broken symmetry. The arrow indi- 

cates a possible choice of vacuum. 



Since QCD describes a very large collection of phenomena at high energies extremely well, there 
must thus be another way to include this symmetry in the real world. This was found by Goldstone [2H] 
and is often called the Nambu-Goldstone mode, while a direct realization is referred to as the Wigner 
or Wigner-Eckart mode. Nambu's papers for this are Ref. |29j . 

Let us first describe this mode for a simpler model. A complex scalar field with Lagrangian 

C = dy*d^ - V{4>) . (22) 

We first look at a potential of the type shown in Fig. ^with a standard form of the type 

V{(j>) = /i 2 0*0 + A (0*0) 2 . (23) 



We choose here A > 
phasetransformation 



to have a stable theory. This Lagrangian has a U(l) symmetry under the 

-► e~ ia (j) . (24) 



This transformation is rotation around the z-axis in Figs. Q and El 

If we choose /z 2 > 0, the potential V(<f>) has the form shown in Fig. ^ where the horizontal axes 
are the real and imaginary part of while the vertical axis are V(<f>). In order to have a full theory 
we have to determine first the vacuum, or lowest energy state, of the system. The contribution of the 
kinetic term, d^cffdy,, is minimized by a constant and spatially homogenous field 0o- From the form of 
the potential, we can see that the total energy is thus minimized for a value of 0o — 0. I.e. (0) = 0. 
Excitations around the vacuum, which give the particle spectrum, have only massive modes with a mass 
m = \x. Things to remark here: The vacuum is unique, i.e. there is only one possible choice of (0). 
There are two massive real modes in the spectrum corresponding to the real and imaginary part of 0. 
The interactions of these particles are simply the four boson vertex directly present in the Lagrangian 
(|2*2*Jl . This mode corresponds to the most standard realization of symmetries like the realization of 
rotation symmetries in standard quantum mechanics. States thus fall in multiplets of the symmetry 
group and amplitudes obey the relations of the Wigner-Eckart theorem. 

However, when we choose the potential with the same form but take /i 2 < the potential looks 
differently as depicted in Fig. El The potential is still invariant under the symmetry (|23j) . but now we 



have more than one option to choose from for the lowest-energy state. The different choices are related 
by a symmetry transformation. In order to start determining the particle spectrum we need to choose 
a particular vacuum state (the lowest total energy is still given by constant due to the kinetic term 
otherwise giving a positive contribution). The arrow in Fig. |21 indicates a possible choice for (0) = 0o. 
All possible choices are of the form 

(0) = 0o = ^f ia » (25) 

with 



-/*7A • (26) 

Different choices of «o lead to the same physics. We can now try to get at the excitations around the 
vacuum. For that we need to parameterize the changes of around 0o- This parameterization can be 
done in many ways but let us choose here the form 

(j)(x) = {v + n{x)) e <(«o+*(*)/«) _ ( 27 ) 

Putting this into the Lagrangian (J22)) we obtain 

£ = \d^d^ + -a 2 n 2 - Xvn 3 - -Xrf - lfi 2 v 2 - -\v A + \d^d^ + \ (2vr) + n 2 ) d^ird^it . (28) 

Z 4r t: Z t: Z Z 

What do we see now, the Lagrangian has no obvious remainder of the symmetry we originally had 
in Eq. (|24|). We say that the symmetry is spontaneously broken. The fact that we needed to make 
a particular choice of the vacuum state means that the original symmetry is no longer visible in the 
spectrum nor in the interactions. There is a massive mode, n, with mass m 2 = — /x 2 /2 and a massless 
mode, it. For the former there is some reminder of the symmetry in the relations of the cubic to the 
quartic coupling and for the latter there is a more striking result. It only interacts in a way that vanishes 
for zero momentum. This is called a low-energy theorem and the massless mode is called the Goldstone 
boson. We also see that the physical results do not depend on the precise choice of the vacuum, a has 
disappeared from the Lagrangian directly describing the excitations around the vacuum of Eq. (J28J1 . 

The appearance of both phenomena, a massless mode and the vanishing of its interactions follow 
from the fact that we have to choose a vacuum in every spacetime point. A qualitative description is as 
follows. The massless mode corresponds to choosing slightly different vacuum states in each space time 
point. This gives a small kinetic energy but also momentum. We can think of this as rolling around in 
the bottom of the valley in Fig. |21 The interaction must vanish at zero momentum since the absolute 
choice of vacuum cannot matter, we can thus shift tt(x) by an arbitrary amount and no physics results 
should change. The interactions can thus only proceed via derivatives, hence the low-energy theorems. 

It is often said that the symmetry is realized nonlinearly in the Nambu- Goldstone mode. This can 
be seen from the fact that the transformation on tc(x) under the original symmetry ({21} is 

ir(x) — > n(x) — iva . (29) 

Note that the new Lagrangian (J2"K|) is invariant under this. We will also talk in the remainder about 
a nonlinearly realized symmetry but will construct for the case relevant for ChPT objects that do 
transform linearly since this simplifies constructing the Lagrangians. 

I have used a very simple first model with the Lagrangian of Eq. fl22(l . to describe the most important 
features here. Let me now shortly indicate which parts generalize to other theories. In general we have 
a continuous symmetry group G which is generated by a number of generators T a . A generic group 
element can be schematically written as g = exp(ie a T a ) . The choice of vacuum leaves in general not 
all generators invariant. The subspace of generators, T b that leave the vacuum invariant generates the 
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Table 1: The two different modes of symmetry realizations compared. 



Wigner-Eckart mode 


Nambu-Goldstone mode 


Symmetry group G 


G spontaneously broken to subgroup H 


Vacuum state unique 


Vacuum state degenerate 


Massive Excitations 


Existence of a massless mode 


States fall in multiplets of G 


States fall in multiplets of H 


Wigner Eckart theorem for G 


Wigner Eckart theorem for H 




Broken part leads to low-energy theorems 


Symmetry linearly realized 


Full Symmetry, G, nonlinear ly realized 




unbroken part, H, linearly realized 



unbroken subgroup H with elements of the form exp(ie b T b ) . The broken part of the symmetry group 
corresponds to those generators that move the vacuum around, i.e., 



T c |0>^0. 



(30) 



The space on which the Goldstone bosons live is the space of possible vacua. This space has the structure 
G/H, the coset space of the group G with its unbroken subgroup H. In our example, this coset space 
was the bottom of the valley and could be simply parameterized by tt(x). The effective Lagrangian after 
spontaneous symmetry breakdown can be made fully invariant also under the full symmetry group but 
this in general implies nonlinear transformations for the Goldstone bosons. 

Alternative introductions to spontaneous symmetry breaking can be found in most books on particle 
physics when the Higgs mechanism is introduced. A review that discusses many more places where 
Goldstone bosons and effective Lagrangian appear is Ref. [3U]- More mathematical descriptions can be 
found in the lectures by Pich jl] and Scherer jSJ |B] . 



2.4 Spontaneous symmetry breaking in the presence of explicit symmetry breaking 

In QCD, the chiral symmetry is not exact, the term with the quark masses is not invariant under chiral 
symmetry. So what happens if there is explicit symmetry breaking present at the same time as the 
spontaneous symmetry breaking. We will again discuss it in the framework of our simple model with 
a U(l) symmetry and a complex scalar field. We now add to the Lagrangian of Eq. (|22jl a term of the 
form 

£extra = ~P(.<P + ^) ■ (31) 



This extra term is not invariant under the symmetry transformation ({24}) and the potential looks tilted 
as shown in Fig. El What is important is that the tilting is small compared to the height of the central 
bump so it can still be treated as a perturbation. The vacuum in Eq. (J25J) is no longer unique but only 
the one with a = is the lowest energy state. 

The physics has also slightly changed. We expand around the vacuum again by using (|2~7j) with 
«o = and obtain in addition to the terms in Eq. l2*%|) 



£ 



extra 



-P>/2(v + rj) (e in(x)/v + e - in{x)/v ) . 



(32) 
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Figure 3: The potential V(0) for a spontaneously broken symmetry in the presence of a 
small explicit symmetry breaking term. The arrow indicates now the only possible choice of 
vacuum. 



The linear term in r\ can be removed by a small additional shift. This happened because the lowest 
energy state is slightly shifted compared to the value v = J— y? j\. But more importantly, when we 
expand the exponentials, we now find that the 7r(x)-field has gotten a small mass, small compared to 
the mass of the 77-field, and no longer has only derivative interactions. The n mass 

9 2^2/3 , s 

v 

is small and can be expanded in the small symmetry breaking parameter (3. The particle corresponding 
to it, is now called a pseudo-Goldstone boson. As long as the explicit symmetry breaking is small, we 
can still use Goldstone's theorem as a first approximation and then add the corrections systematically. 
This is precisely what we do in ChPT when the light quark masses are explicitly included. 

2.5 Spontaneous symmetry breaking in QCD 

We already argued in Sect. 12.31 that the chiral symmetry of QCD cannot be realized in nature since 
the predicted parity doublets do not occur. We thus expect the chiral symmetry to be realized in the 
Nambu-Goldstone mode. What theoretical evidence do we have directly for this? 

Most of the remainder of this paper is about the Goldstone bosons from the spontaneous chiral 
symmetry breakdown and their properties. In this way, all those properties are strong indications that 
the picture described below is correct. However let us first give the full theoretical arguments. 

• It has been proven that the chiral symmetry is spontaneously broken in the limit of a large number 
of colours and assuming confinement [SI] . 

• The vector symmetries remain unbroken in a vectorlike symmetry as QCD |32j . 

• Assuming confinement, the anomalies in the effective low-energy theory must match those for the 
underlying QCD theory. For two flavours, this can be done but not for three or more flavours. 
We thus need spontaneous symmetry breaking in order to have a correct anomaly matching for 
three or more flavours [33J. 

We thus believe that the flavour symmetry SU(np) x SU(np) is spontaneously broken down to the 
diagonal subgroup SU(tif)v — SU(tif)l+r also for the realistic case of three flavours. There are eight 
broken generators and we thus expect eight Goldstone boson degrees of freedom. If we look at the 
hadron spectrum there are eight natural candidates for this. The three pions, 7r°, ■k ± , four kaons, K ± , 
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K°, K° and the eta, rj. Goldstone's theorem thus explains their low mass compared to the other hadrons 
as well as the fact that their interactions are relatively weak. The three main qualitative predictions 
that follow are: 

1. The masses are rather small. 

2. The masses obey the Gell-Mann-Okubo relation with squared masses rather than linearly. 

3. 7T7T scattering is fairly small compared to proton-proton scattering and related to pion decay. 

4. There is a nontrivial relation between the pion decay constant, the axial-vector coupling of the 
nucleon g^ and the pion nucleon coupling G n , the Goldberger-Treiman relation J33]. 

All of these predictions are well borne out by experiment. The first three will be discussed in detail 
below. 

In Sect. l2.3l the U(l) symmetry was broken by a vacuum expectation value of the field <fi. In QCD, the 
vacuum is also not invariant under the full chiral symmetry G x but the quantity that most characterizes 
the noninvariance of the vacuum is a composite of the fundamental fields, the quark-antiquark bilinear 
condensate. The vacuum of QCD is thus characterized by 

(^)^O. (34) 

This is the standard picture which we now know is true for the two-flavour case jSH]- Everything we 
know indicates that it is also true for the three-flavour case but the argument is not fully closed yet [36J . 
That the vacuum expectation value (|34j) breaks chiral symmetry can be easily seen when we rewrite 
it into left and right handed components and as a matrix in the flavour space via 

(V% = (q jL q lR ) . (35) 

Under a chiral symmetry transformation g^ x g R e SU{np)p x SU(np) R this transforms as 

V q - g R V q g{ . (36) 

We now choose a particular vacuum expectation value, 



(V% = -5 rj (qq) , (37) 



where we used the generic symbol (qq) for a vacuum expectation value of one quark species in the 
chiral limit. With this choice, one sees that the vector subgroup, g^ = g R , leaves V q invariant but 
any transformation with gi ^ g R does not, so the axial part of the symmetry group is spontaneously 
broken. 

It now remains to parameterize the space of possible vacua. For a spontaneous breakdown of 
G = SU{uf)l x SU{uf)r to H = SU{uf)v the coset space G/H has itself the structure of an SU(rip) 
manifold. The simplest parameterization for describing the Goldstone boson is thus choosing them as 
an rip x np special unitary matrix U . This parameterization, together with some minor extensions, is 
used in the remainder of this paper. 

We argued in Sect. 12.31 that in some sense the Goldstone bosons live in the space of possible vacua. 
The same is true here. We can parameterize the space of vacua of V q by the same special unitary 
matrix U via 

V q = \u(qq) ■ (38) 
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The matrix U thus transforms under the chiral symmetry group as 

U^g R Ug L . (39) 

There exist many possible alternative parameterizations. The solution for the two-flavour case was 
originally found by Weinberg [SZj. It was generalized by Coleman, Wess and Zumino to arbitrary 
symmetry breaking patterns G — > H 38J. The inclusions of states other than the Goldstone bosons 
was worked out in Ref. [HH]- The latter shows also explicitly that there is no need for the existence of 
states related by parity in order to have a fully chirally symmetric theory when the chiral symmetry is 
spontaneously broken. The last two references also showed that their parameterization is fully general 
and remains valid when loop effects are taken into account. 

2.6 Lowest Order ChPT 

Chiral Perturbation Theory is the low-energy effective field theory of QCD where the degrees of freedom 
taken into account are the Goldstone bosons from the spontaneous breakdown of the chiral symmetry 
and their interactions. We showed in Sect. 12.51 that the resulting Goldstone boson manifold can be 
parameterized by a special unitary matrix U which transforms under the chiral symmetry group as in 
Eq. ()39|). We also want to include the external fields introduced in Sect. 12.21 and we want it to be fully 
invariant under the chiral symmetry as required by the arguments of [3KJ- Note that these arguments 
including the loop level are worked out in great detail in Ref. |4T)] . 

In the remainder a lot of notation will be introduced. In particular many traces of np x np matrices 
will appear. In order to make these traces easier to see we introduce the notation 

(A) = tip (A) , (40) 

where trp denotes the trace over flavour indices. 

Without external fields and derivatives the only terms that can be constructed are of the form 

C = a (U j U}+a 1 detU + a* l detUK (41) 

This is only an irrelevant constant since for a special unitary matrix we have that 

U ] U =1 and det U = 1 . (42) 

This corresponds to the fact that Goldstone bosons cannot have interactions without derivatives or 
explicit symmetry breaking. 

At the next order, we find that d^U is not chirally invariant. A building block that transforms nicely 
can be constructed by defining a covariant derivative 

D^U = dpU-irnU + iUlu, 
D,JP = d^-il^ + iU^r^. (43) 

These transform simply under the chiral symmetry group as 

D,U - g R D^Ug{ Djfl - g L D^g R . (44) 

With the covariant derivatives in hand we can now construct the Lagrangian at the first nontrivial 
order 

C 2 = /3 1 (D lt rtD' t U) + f3 2 (MU ] + UM ] ) + i(3^{MU ] - UM ] ) . (45) 

In this construction we have assumed that 

(r>> = (U = . (46) 
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This together with det U = 1 implies that 

(rfD (t U) = 0. (47) 

We have used the fact that the Lagrangian can be changed by partial integration and that partial 
integration can be done with the covariant derivatives. 

Under parity we interchange left and right. The parity transformation on U is thus U — > U* . A4 
goes similarly to its complex conjugate. The term proportional to (3% thus violates parity and can be 
dropped because of that. 

The Lagrangian (|45J1 is usually written in the form 



F 

with 



C 2 = ^-(D^D^U + X U ] + U X ] ) , (4* 



X = IBM = 2B (s + ip) . (49) 

The specific values of the parameters F and B depend on the number of flavours and they are conven- 
tionally written as F, B for the two-flavour case |2] and F , B for the three-flavour case [3]- 

The Lagrangian for lowest order ChPT was first given by Weinberg in Ref. |37| for the two-flavour 
case and soon generalized. A discussion about its construction, including many of the subtleties in- 
volving the overall phase of U can be found in Sects. 3 and 4 in Ref. [3] and the references mentioned 
therein. 

2.7 A few consequences of lowest order ChPT 

In this section I will only talk about the three-flavour case, rip = 3, and thus stick to the notation for 
that case. Let us derive a few simple consequences from the Lagrangian (|48j). This can be done using 
the machinery of the effective action or by simply using Feynman diagram calculations. Both methods 
have to give the same answers and calculations have been performed using both approaches. In this 
section we will stick to the simplest one, tree level Feynman diagrams and identifying masses by the 
terms in the Lagrangian. 

First we need to parameterize the matrix U. This is done in terms of a matrix of meson fields in 
the form 

(JS^ + JSV -+ K+ \ 



U = e iV2M/F with M 






(50) 



V A AU -75 V 

We have used here the isospin triplet field for the 7r° and the octet component only for the r]. This is 
fine as long as we work in the isospin limit with 

m u = m d . (51) 

In this limit we have two quark masses 

rh = - (m u + m d ) (52) 

and m s . We use the relation (jI7j) and set s in (|4"8|) equal to M.. We put the parameterization of U into 
(I48JI and expand the exponentials. Looking at the terms up to second order in the meson fields, we find 

£ 2 = -<9 M 7T°<^7r + d^K+d^K- + d^d"K° + d^n+d^n- + -d^d^ - B m n\° - 2B m u+ij- 

-B (m + m 8 ) K°K° - B (m + m s ) K + K~ - B —- — - r)r) . (53) 
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Figure 4: The Feynman diagram responsible for the main pion decay n + — > fi + u. 

Here we see several things. The pions have the same mass 

m\ = 2B m . (54) 

Similarly, the kaons have the same mass, 

m\ = B (m + m s ) . (55) 

A relation between the pion, eta and kaon masses exists. This relation is the famous Gell-Mann-Okubo 
(GMO) relation: 

m l = ^ m K ~ ^ m l ■ (56) 

We see here that this relation should be satisfied by the masses squared. A naive application of the 
Wigner-Eckart theorem and the symmetry group SU(3)v would have led to the same relation but 
with the masses present linearly. The fact that it, t] and K are pseudo-Goldstone bosons from the 
spontaneously broken chiral symmetry explains why the relation should be with the quadratic masses. 
This follows if we include the additional assumption that the lowest order term gives the bulk of the 
observed masses, see e.g. Ref. |4*Tj . 

The determination (J51J1 of the pion mass actually contains another famous relation. We can get the 
quark- ant iquark bilinear condensate by taking functional derivatives of the generating functional G. 

(uu) = - (-^t-.g) , (57) 

evaluated at the point where the external fields are set to zero. Doing this we obtain to lowest order 

(qq) = {uu) = (dd) = (ss) = -B F$ . (58) 

We thus obtain the celebrated Oakes-Renner relation 



2 m(uu + dd) 

w 



mi = ^ '- . (59) 

o 



A third example is the semileptonic decay of the pion. This proceeds via the diagram shown in 
Fig. EJ The ii + W~ vertex can be derived from the Lagrangian ()48j) when we use the way to include 
the W-boson via (j21j) . Alternatively, the coupling of the pion to the iy-boson is regulated by the pion 
decay constant defined by 

(0|d7 M 7 5 u|7r + (p)) =i^ Pfl . (60) 

We can compare the two calculations or calculate the matrix-element (160(1 directly by taking a functional 
derivative with respect to a^. In both cases we reach the lowest order result 

F„ = F . (61) 
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Figure 5: The Feynman diagram responsible for pion-pion scattering at lowest order. The 
lines are the pions. 



The final result we will allude to is mr scattering. This was first derived by Weinberg using current 
algebra methods J32] • Here it follows from expanding the Lagrangian (|48|) to higher orders in the meson 
fields. We then find vertices depicted schematically in Fig. |S| In terms of the amplitude A(s,t,u) 
defined later in (jl03|) the result he obtained is 

A(s,t,u) = -^(s-ml) . (62) 

■K 

2.8 Powercounting and renormalization overview 

The main purpose of this paper is to review higher order calculations in ChPT. The Lagrangian in 
(|4*8Jl is nonrenormalizable. This makes it unfit to be used as a fundamental theory but produces no 
problems for effective theories. We know that there is a well-defined gauge theory, QCD, underlying 
ChPT. But still, in order to have a phenomenological usefulness, there must be a way to limit the 
number of parameters that are present. In general, nonrenormalizable theories have an infinite number 
of adjustable parameters. 

What we will show in this subsection is that there exists a well-defined way to order the various 
contributions of in terms of expansion parameters. First, there are many quantities here, and we are not 
simply performing an expansion in a small coupling constant as is done in Quantum Electrodynamics. 
The expansion we have here is a long-distance or a small momentum expansion, together with an 
expansion in the quark masses. We first call the magnitude of a typical momentum component p. Since 
these come from derivatives, it is natural to also take the external fields l^ and r^ of order p, since 
these occur together with the derivative in the covariant derivative of Eq. ()43j) . On-shell particles have 
p 2 = m 2 . It is therefore natural from ()54j) to take a quark-mass as order p 2 . This in turn makes it 
natural to count scalar and pseudoscalar external fields as order p 2 because of the rule used in (jT7j) . 

With this counting we see that all terms in the Lagrangian (|48|) are of order p 2 , which is why we 
chose a subscript 2 there. This counting can be generalized to all orders and is how we will order our 
series. It was introduced by Weinberg in Ref. pQ. 

Let us first give it for a few simpler diagrams. On the left hand-side in Fig. |Hl the rules of counting 
are shown. A vertex from the lowest order Lagrangian ([48)) counts as order p 2 . A propagator is of 
order 1/p 2 and a loop integral is of order p 4 . For dimensional reasons it must give an extra four powers 
of momenta. On the right-hand side we show two loop diagrams contributing to pion-pion scattering. 
They are both order p 4 when counting the number of loops, vertices and propagators. The two one-loop 
diagrams are thus of the same order and also of the same order as a tree level diagram with a vertex 
with four derivatives would be. 

This type of observation is the underpinning of the expansion used in ChPT. Let us know give this 
argument in general. A generic diagram with Np propagators, Nl loop integrations and N n the number 
of vertices of order p n . The total order of the diagram is then 

N T = -2N P + 4N L + Y, nN n . (63) 



n>2 
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Some diagrams 
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(p 2 y(i/p 2 ) 2 P 4 = p 
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fd*p p 4 (yr)(l//r )/ = />' 



Figure 6: The power-counting introduced by Weinberg illustrated on the example of pion- 
pion scattering. See text for explanations. 



Here we already used the fact that the lowest order Lagrangian is of order p 2 . We can rewrite this using 
the relation between the number of internal lines, Nj = Np, the number of vertices Ny and the number 
of loops N L , 

Nj = N L + N v -l. (64) 

Since Ny = J2 n N n , Eq. (}63|) can be rewritten as 

N T = 2 + 2N L + J2(n-2)N n . (65) 

n>2 

Eq. (|65|) is the basis of the perturbative expansion of ChPT. The lowest order contribution to any 
process is given by a tree level diagram with only vertices from C 2 - The next order, NLO, is formed 
by one-loop diagrams with only vertices from C 2 and tree level diagrams with vertices from £ 2 and 
one-vertex from the p 4 Lagrangian £ 4 . 

Eq. (J63j) also shows the importance of Goldstone's theorem for the existence of a perturbative 
expansion. Only because the lowest order has derivatives or external fields do we have an expansion 
where higher loops imply higher powers of p, i.e. Goldstone's theorem is the source of the requirement 
n > 2. 

Note that the powercounting described here is closely related to the notion of superficial degree of 
divergence described in most Quantum Field Theory books see e.g. |2o1 I43j . 

The relation given in Eq. (jMj) can most easily be understood by induction. In a tree level diagrams 
all vertices need to be connected by internal lines. The simplest diagram has one vertex and no internal 
line. Keeping it at tree level but adding lines implies always adding one internal line and one vertex, so 
far tree level diagrams we have Nj = Ny — 1. Every time we add an internal line, but no new vertex, 
we create a new loop. We thus end up with (JMJ) . 

Let me finish this section by giving the overview and general arguments involved in ChPT and its 
construction and renormalization. First we use Weinberg's conjecture [Tj, see also the discussion in |44j : 
if one writes down the most general possible Lagrangian, including all terms consistent with assumed 
symmetry principles, and then calculates matrix elements with this Lagrangian to any given order of 
perturbation theory, the result will simply be the most general possible S-matrix element consistent with 
analyticity, perturbative unitarity, cluster decomposition and the assumed symmetry principles. 

Then we assume that the relevant degrees of freedom are the Goldstone bosons from the spontaneous 
breaking of chiral symmetry and construct the most general Lagrangian with them which has the full 
chiral invariance. Using the results of Ref. [3H] this can then be brought into a standard form. We have 
assumed here that we can use a local Lagrangian. That this can be done was shown in Ref. jJO] • As 
said above, Goldstone's theorem implies that the lowest order Lagrangian is of order p 2 . We now use 
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a regularization that conserves chiral symmetry. Dimensional regularization [12] is the most standard 
choice. In a general Quantum Field Theory with local vertices, all divergences that appear are local. 
Since we start with a Lagrangian invariant under the symmetry and a fully invariant regularization, 
the divergences are local but will have a structure that obeys the symmetry structure. Since our 
constructed Lagrangian includes all possible terms consistent with the symmetry, all divergences can 
thus be absorbed into the coefficients of the Lagrangian. The total number of local terms in the 
Lagrangian will be infinite, but since we can order the expansion in terms of the order in powercounting 
in p, we have a well-defined system with a finite number of parameters up to any given order in p. 
A much more extensive version of this discussion where much attention is paid to all the issues just 
mentioned here is Ref. [30] . 

It is possible to use a regularization which is not chirally invariant. One then needs to introduce also 
non-invariant counterterms and explicitly enforce all Ward identities. Some of the problems involved 
are discussed in the papers listed in Ref. [46]. 

Renormalization is a wide topic and can be found treated in most Quantum Field Theory books, 
e.g. [123 ESI, one concentrating on renormalization is |47j . 

2.9 Construction of higher- order Lagrangians 

In Eq. (J48J) we showed the lowest-order Lagrangian. In Sect. 12.81 we presented how ChPT can be 
systematically extended to higher orders. A major part of this involves constructing the most general 
Lagrangian at a given order in the powercounting in p. This involves two steps. First we want to 
construct a complete Lagrangian that includes all possible local terms that are invariant under the full 
chiral symmetry. This is a rather elaborate exercise but can be done in a fairly straightforward manner. 
However, this procedure tends to end up with far too many terms. The more challenging part is to find 
a minimal but still complete set of terms. 

To construct a complete Lagrangian one first constructs a complete set of quantities involving U, 
derivatives and the external fields that transforms in a simple manner under the chiral symmetry. As 
an example, for a quantity with one derivative there are three standard choices 

L^ = iU ] DJJ R^ = iUD^U ] and u^ . (66) 

The quantity u^ needs a little more explanation. We write the full matrix 

U = u 2 . (67) 

For a general chiral symmetry transformation gi x g^ g SU{np)L x SU{np)R there exists a matrix 
h G SU{n F ) v such that 

u -> g R uh) = hug{ . (68) 

Eq. ()68|) is the definition of h. The unitary matrix h depends nonlinearly on g^ g R and u , but is 
unique. This is really the general parameterization of [3HHHH] for the case of SU(n) x SU(n) — > SU(n). 
We now define 

u^ = i [u\d^ - ir M ) u-u(d^- il^) u f j . (69) 

Under the symmetry, the three choices transform as 

L^ —> gL^^L j 
Rfj, —> 9rR^9r , 
?V ~~> hu^h) , (70) 

which follow from Eqs. (J39I44I15I70J) . The kinetic term of the lowest order Lagrangian can be written 
in terms of all three using 

(L^) = (R^) = {uX) ■ (71) 



T M is the connection 



In this review we use the last choice but the others have also been used, see e.g. ^SJ EH EDI- 
To get the order p 4 Lagrangian, we need the additional quantities 

X± = u ] xu ] ±ux ] u, 

/I" = uFr^±u^F^u, (72) 

where F L and F R denote the field strengths of the external fields / and r, such that 

Ff = 3T - d v F - i [F, l v \ . (73) 

F^ is defined analogously in terms of r. All the quantities in Eq. (|72*|) transform under chiral symmetry 
as 

X -> hXh) . (74) 

In f!43p we defined a covariant derivative that transforms simply. For objects transforming as (|74*|) a 
covariant derivative can also be defined via 

V fl X = d fl X + T ll X-XT fl . (75) 

r M = - {w f (<9 M - ir^) m + m (^ - Up) m 1 "} . (76) 

Using the transformations defined earlier, it can be shown that V M A" transforms as (ff4*|) as well. One 
last relation which can be checked by putting in all definitions is: 

f-tiv = V„tt M - V^Uy . (77) 

The most general Lagrangian of order p 4 after using partial integrations and all the identities men- 
tioned above for the case of rip flavours is OH] 

12 

£4 = J^ LiXi + contact terms 

= L (u^u u u„u v ) + Li (w%) 2 + L 2 (uV) (u^u v ) + L 3 ((w%) 2 ) + L 4 (m\) (x+) 

+L 5 («%*+} + U (X + ) 2 + L 7 (x-) 2 + y (xl + X 2 -) - %U (/TV**) + *f (f + - f_) 
+ii*i (x- (V% - lx-)) + L12 / (V% - ^-) 2 \ + H x (Fl + F 2 R ) + H 2 (jo*), (78) 

where the definition x_ = X- — (X-)/ n F has been applied. Furthermore, the lowest order equation of 
motion is given by 

X EOM = V% - l -X- = 0. (79) 

The Lagrangian of Eq. (|78j) contains three types of terms. The last type, the terms proportional 
to Hi are contact terms which contain external fields only. Thus they are not relevant for low-energy 
phenomenology, but they are necessary for the computation of operator expectation values. Their values 
are determined by the precise definition used for the QCD currents, and they are conventionally labeled 
hi and H\ for unquenched %PT with nj = 2 and n/ = 3 quark flavors, respectively. 

The terms containing L\\ and L\i are proportional to the equations of motion. They can thus always 
always be reabsorbed into the Lagrangians of higher orders. A full proof can be found in Ref. |51j . 
App. A. A discussion at lower level is Ref. [S2J. The fact that these can be put into the Lagrangians 
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of higher order by field redefinitions can be qualitatively understood by the following argument. The 
variation of the lowest order Lagrangian gives the equation of motion. Using U = ue^u for the variation 
of U and varying £ gives 

5C 2 ex (£X EOM > • (80) 

A term in a Lagrangian of the form (AXeom) can thus removed by a field redefinition of the type 
U — > ue %A u . The derivation with all factors correct and worked out to all orders can be found in 
App. A of Ref. [SI] . 

The physically relevant terms are the remaining ones, containing L^ i — 0, . . . , 10. Just as for F 
and B, the Lj are different for every value of the number of light flavours rip. 

For a specific number of flavours, additional relations exist, the Cayley-Hamilton relations. These 
follow from the fact that any n- dimensional matrix satisfies its own characteristic equation. This is 
described in Sect. 3 of Ref. [5T] . For the p 4 Lagrangian for three flavours this allows for the removal of 
the term proportional to L via 

(uVvO = -2<u%u"m i/ ) + -(wX)(uX> + (uVX^u.) . (81) 

This is the same as Eq. (7.24) in Ref. [Sj. The Cayley-Hamilton for two flavours is 

{A,B} = A(B) + B(A) + (AB)-(A)(B) (82) 

for arbitrary 2x2 matrices A, B. An additional contact term exists for two-flavours as well at order p 4 
since det x is invariant under chiral SU{np) x SU(tif) transformations and is order p 4 for two flavours. 
As a result there are 7 U and 3 hi parameters at order p 4 for he two-flavour case [2] . The correspondence 
with the general number of flavours Lagrangian is via 



h -- 


= -2L + 4L 1 + 2L 3 , 


k -- 


= 4L + 4L 2 , 


k -- 


= -8L 4 -4L 5 + 16L 6 + 8L 8 


k -- 


= 8L 4 + 4L 5 , 


h -- 


- L 10 , 


k -- 


= -2 Lb, 


h -- 


= -16L 7 -8L 8 . 



(83) 

We can for the two-flavour case write the Lagrangian in the form (|78|) but only the combinations in 
(JSHj) will show up in experimentally relevant quantities. 

All the same principles apply to the construction of the Lagrangian at order p 6 . However, since many 
more combinations are possible, it becomes much harder to find a minimal set. This was accomplished 
in Ref. J3TJ after the first attempt of Ref. J2HI- We will not show the full Lagrangian here, it is given 
in the appendices of Ref. [SI]. There also all the Cayley-Hamilton relations that were used to obtain 
the minimal set for the two- and three-flavour case can be found. In Tab. 01 the number of independent 
parameters at each order in the Lagrangian is summarized for the cases relevant in this review. 

The Lagrangians at order p 6 contain very many terms and the relevant operators can be found in 
Refs. |51| HkI] . The parameters are labeled Ki for the general-ni? flavour case, Cj for the three-flavour 
and Cj for the two-flavour case. 

The parameters in the Lagrangians are often referred to as low-energy constants (LECs) and the 
terms in the Lagrangians are sometimes called counterterms even though the latter strictly only means 
the additional divergent parts defined in Sect. 12.101 . We will use the terms parameters in the Lagrangian 
and LECs interchangeably. 

20 



Table 2: The relevant sets of LECs, where the i+j notation denotes the number of physically 
relevant (i) and contact (J) terms in the respective Lagrangians. At NLO, the latter ones 
are conventionally denoted hi for n/ = 2 and Hi for n/ = 3. 



n f 


XPT 
2 


XPT 
3 


XPT 

n 


PQxPT 

2 


PQxPT 

3 


LO 


F,B 


Fo,B 


F ,B 


F,B 


Fq,Bq 


NLO 

i+j 


7 + 3 


10 + 2 


Li 

11 + 2 


L (2pq) 

11 + 2 


iJZpq) 
11 + 2 


NNLO 
i + j 


Ci 

53 + 4 


Ci 
90 + 4 


Ki 

112 + 3 


K (2pq) 

112 + 3 


js(3pq) 
112 + 3 



2.10 Renormalization in practice 

In this review we will not treat renormalization in detail. A short overview is given in Sect. 12.81 A more 
comprehensive discussion, including more details, of renormalization in ChPT can be found in Ref. j54j . 
A general treatment of renormalization is the book by Collins jUj • 

In ChPT one uses in general dimensional regularization to regularize the divergences that occur. 
The number of space time dimensions becomes noninteger and is written as 

d = 4 - 2e . (84) 

All integrals are expanded in a Laurent-series in e and the divergences occur as inverse powers of e |33] . 
These so-called pole-terms are then absorbed by adding to the parameters in the Lagrangian also terms 
divergent when e is sent to zero in such a way that the final result is finite. When one only adds terms 
sufficient to precisely cancel the poles the procedure is called minimal subtraction (MS). But, there 
are several extra pieces that always show up together with the poles [53]. These can be subtracted 
together by adding finite parts into the additional terms 1^3, a procedure known as modified minimal 
subtraction (MS), corresponding to choosing c 7^ 1 in Eq. (}85j) below. 
To order p A at most single poles can appear and we define 



Uid) 



{tic) 



-2c 



(47T) 



~ + L' i {fi,c,e) 



The usual choice in ChPT is J2] 



1 



lnc=--[ln47r + r'(l) + l] . 



35) 



(86) 



Note that the renormalization procedure in Eq. (|85j) has also introduced a scale //. The usual choice of 
low-energy constants is 

L r M = Ll^c,0). (87) 

One more subtlety is the choice of the order e part in 

L r i (ji,c,e) = L r i ( f M 1 c,0)+6 i e + ..., (88) 

in this review we choose Si = 0. Different choices are possible but the difference can always be reabsorbed 
in the values of the parameters at higher orders 
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At order p 6 double poles can occur and the equivalent of (|85|l reads |54*1 lo^] 
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The coefficients f j are known for the 2,3 and iV^-flavour case El EH] arid likewise the coefficients 
I] , Tj and T\ [S3]. These coefficients can be calculated in general using heat-kernel techniques and 



background field methods, see Ref. [56] for the methods and earlier references. 

Note that the renormalized K[ have been made dimensionless by the extra factors of F, while the 
Ki have dimensions. This extra factor of F 2 has not always been treated consistently in the literature. 

The renormalized coefficients obey renormalization group equations and there are consistency con- 
ditions between the various coefficients appearing in (|85j) and (|59*|) . These were first discussed by 
Weinberg [T] and are called Weinberg consistency conditions. They were used in Ref. [23 EH] to obtain 
first two-loop results. In the notation introduced above they are 
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The last equation in (|9"T|) implies that the coefficient of the double pole can be derived from only one-loop 
diagrams PJJ EZ1 EH] • This has been extended to higher orders in Ref. [SHj • 

The renormalization coefficients Tj needed in (|85jl for the cases of 2,3 and uf flavours are denoted 
by 7j, Tj and fV The two-flavour ones are [2] 



7i = o5 ^2 = -; 73 



The three-flavour ones are [3] 
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while the rip flavour case is given by |3J 
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Figure 7: The imaginary part of the pion-form-factor to two-loop order. It only needs the 
pion-pion scattering and the pion form-factor to one-loop. The lines are pions, the double 
line is the insertion of the external field. The thick dots indicate a general diagram, the 
small dot a vertex from £2- Adapted from Ref. jHO] 



For the two-flavour case numerical values are usually not quoted for l\ at a given scale but instead 
in terms of the barred quantities defined by j2] 
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These are independent of the scale \x. 



3 Two-flavour calculations 

3.1 Using dispersive methods for p 6 

The 5-matrix is unitary. The T-matrix defined by S = 1 + iT, thus satisfies the relation 

T ] T = -i(T - T f ) . 



(96) 



This leads to various relations between the real and imaginary part of amplitudes. These can be brought 
into relations for each diagram separately and several parts of diagrams can thus be constructed from 
their imaginary parts. The latter are well defined via Cutkosky rules. The first phenomenological 
applications of ChPT beyond order p 4 were done using methods of this type. In Ref. [HE] this was used 
to determine parts of the p 6 corrections to pion form-factors. The parts that can be determined this way 
are the nonanalytic dependences on kinematical variables. With kinematical variables here is meant the 
momentum transfer q 2 for form-factors, the Mandelstam variables s, t, u for scattering processes and 
their equivalents for other processes. 

The method works by evaluating the imaginary parts using the Cutkosky cutting rules. The real 
parts can then be worked out using Cauchy's rule when sufficient subtractions are made to make the 
dispersive integral convergent. Due to the subtractions the analytical dependence on the kinematical 
variables cannot be reconstructed in this way. The principle is illustrated in Fig. [7| We can see there 
that the knowledge of pion-pion scattering and the form-factor to one-loop is sufficient to determine 
the imaginary part of the form-factor to two-loop order. 

The methods used in JH01 have been extended to other processes as well. In particular, the structure 
of ChPT together with the Roy equations was used to calculate the kinematical dependences of pion- 
pion scattering to order p 6 in Ref. [SI] . The observation made there, which has since been generalized to 
many similar processes, is that the imaginary part up to order p 6 in scattering processes only depends 
on at most one kinematical variable in a nonanalytic fashion at a time. This means that up to order 
p 6 the amplitudes can be written in terms of single-variable functions. This has been a very useful 
observation in simplifying many of the full order p 6 calculations done afterwards. The same reference 
also observed that all needed integrals for pion-pion scattering could be performed analytically. 

The third set of processes to be discussed fully in this way has been the decay r —> 2tc, 37r in Ref. 
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Figure 8: The cross-section for 77 — > ir 7r° with | cos#| < 0.8. The data points are from |HH] . 
The dashed line is the order p 4 result from [BUI IM] . The full line is the order p 6 calculation 
from jHZ]. The dash-dotted lines indicate the band of results from the dispersive results [66 . 
The improvement of ChPT in going from p 4 to p 6 is remarkable. Figure from Ref. [6*7] . 



3.2 The process 77 



ttV 



The first process to be fully calculated to two-loops was in fact a rather difficult one. It was the process 
77 — > ir 7T°. The interest in this process for ChPT started because it was realized early on that, to 
order p 4 , this process did not depend on any of the order p 4 LECs JH3J E3|- It thus provided a clean 
prediction from ChPT. This process was then measured by the Crystal Ball Collaboration |J3]. The 
overall size of the prediction of Refs. [63, 64 j was in good agreement with the data but the rise with 
center-of-mass energy predicted by order p 4 ChPT was not seen in the data. 

This fact was repeatedly used to emphasize the inadequacy of ChPT, see e.g. the discussion in 
Ref. jHSl- This prompted the authors of Ref. [U7| to start the calculation of this process to order p 6 . 
It turned out that many of the relevant integrals were in fact not known despite many years of two- 
loop calculations in other circumstances. The necessary techniques were developed by the authors of 
Ref. [68J. The full order p 6 calculation gave a significant improvement over the p 4 calculation as is shown 
in Fig. |H1 taken from Ref. [H7J. It should also be noted that the convergence of ChPT is reasonable in 
the entire range below 600 MeV as can be seen from the figure. 

This calculation has since been redone with improved techniques for the integrals [HH] where the 
results of the earlier calculation have been essentially confirmed. 



3.3 The pion mass and decay constant 

The simpler observables, the pion mass and the decay constant were calculated somewhat later. First 
in Ref. [701 [7T] and later confirmed by [72133]. For the equal mass case, all relevant integrals can be 
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done explicitly. We quote here the result rewritten in terms of the physical mass [73*] . 
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The constants r r F and r r M denote the contributions from the 0(p 6 ) Lagrangian after modified minimal 
subtraction and are given by 



r r F = 8c r 7 + 16^ + 8cg 

r r M = -32c r e - lQc r 7 - 32c r 8 - lQc r Q 4 l.Sr' l0 4 Wc' u + 32r', 7 - G U- 

In Eqs. (j9*7j) and (|55jl we have used the quantities 
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M 2 being the lowest order pion mass and F the pion decay constant in the chiral limit. The l\ are 
the finite part of the coupling constants U in £4 after the MS subtraction as given in Eq. (|85|). The 
ki explicitly include the relations between double logarithms and single logarithms that follows from 
Weinberg's consistency conditions and were first introduced in [q'H] . 



3.4 The process 77 — ► tt + tt and polarizabilities 



The process 77 — > tt + tt~ has been the focus of a large amount of theoretical and experimental attention 
as well. The order p 4 expression was worked out in Ref. [HB! and lead to a clean prediction for the 
polarizabilities, see e.g. |71j . 

The agreement with the neutral pion cross-section and polarizabilities is reasonable as discussed 
in [Un] and [73], especially when the order p 6 corrections corrections are taken into account. References 
to earlier theoretical and experimental results on the neutral pion polarizabilities can be found in those 
two papers. 

For the charged pion polarizabilities, the comparison with experiment is not that good. The order 
p 6 corrections were calculated in Refs. [SHE] and were of the expected size. With the new experiment 
[To] and dispersive derivations, Ref. J7.7J and references therein, there is a significant disagreement with 
the predictions of ChPT. The order p 6 calculation has been redone recently J7H] and is essentially in 
agreement with the older calculation as well. 
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Figure 9: The cross-section of 77 — > tt + tt~ with | cos#| < 0.6 in the regime where ChPT is 
fully valid. The data are from Mark-II +citedataggpipi, the lowest order result is equivalent 
to scalar QED, the order p 4 result and order p 6 are from [H2] and [7TJ [TKj. Figure from 
Ref. US]. 



The agreement with the measured cross-section data is quite good and the ChPT series shows good 
convergence going from order p 2 to p 4 and p 6 . This is shown in Fig.[§]where the theoretical cross-section 
is compared with the data of Ref. [75] • 

The values of the polarizabilities are not in good agreement with the ChPT predictions, but there 
is a wide spread in the direct experimental values and the dispersive estimates. A discussion with more 
references can be found in J78J and [77] . Here I only quote the final result of order p & ChPT [75] 



(«i - AU = (5.7 ± 1.0) • HT 4 fm 3 



and the result from the latest experiment [TE] 



(ai - &U = (11.6 ± 1.5 stat ± 3.0 syst ± 0.5 mod ) • 10~ 4 fm 3 . 



(101) 



(102) 



One sees that there is a clear, at present not understood, discrepancy. The ChPT prediction is not 
expected to have large corrections, it converged well and has good agreement with the direct data on 
77 — > 7r + 7r~. The experiment of [7H] has also large contributions from the direct process jN — > 7./W 
and there might be ununderstood effects in the separation from the pion pole contribution. 



3.5 Pion-pion scattering 

Pion-pion scattering has received an enormous amount of attention both on the theoretical and exper- 
imental front. It is in some sense the most pristine of hadronic processes, involving only the lightest 
strongly interacting state. As mentioned earlier, the calculation by Weinberg in current algebra was one 
of the great successes of that approach in explaining the relative smallness compared to other strong 
interaction cross-sections. 

The amplitude for pion-pion scattering can be written as 

(AP^ C (P3) outKGoOvrV) in) = (n d ( Pi )n c (p 3 ) in|7r a ( Pl )7r b (p 2 ) in) 



+i(2ir) A 5 4 (P f - Pi) {5 ab 5 cd A(s, t, u) + cycl.} (103) 
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where s, t, u are the usual Mandelstam variables, expressed in units of the physical pion mass squared 

M 2 



s = ( Pl + P2 y/M^ t = ( P3 - Pl y/Mt, U = (p 4 - Pl y/M^ 



(104) 



Using these dimensionless quantities, the momentum expansion of the amplitude amounts to a Taylor 

series in 

M 2 
** = -p£> (105) 

IT 

where F n denotes the physical pion decay constant. 

The lowest order result was found by Weinberg using current algebra methods J221) the order p A 
calculation was performed by Gasser and Leutwyler [HU] • The full calculation was done in [721 03] • The 
result obtained there is 
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Figure 10: The pion-pion phase-shift difference <5{J — 5\ as a function of the center-of-mass 
energy. Shown are the lowest order, order p A and order p 6 . The latter is shown for two sets 
of input parameters together with the then most precise data jHT]. Figure from Ref. |54j . 



The loop functions J and Ki are 
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The functions s~ x J and s~ x Ki are analytic in the complex s-plane (cut along the positive real axis for 
s > 4), and they vanish as \s\ tends to infinity. Their real and imaginary parts are continuous at s = 4. 
The coefficients bi in the polynomial part are given in App. D of [S^ and their dependence on the p 6 
LECs can be found in [o^j . 

A few comments are in order here. The result (jl06|) can be rewritten in the form derived in jHJ. The 
relevant loop-integrals can all be written in terms of elementary functions. This is a feature of many 
of the results when all masses are equal but not for all. The ChPT series in fact converges reasonably 
well as can be seen in Fig. E3 

This is a review on ChPT but some of the recent history in the theoretical treatment of pion-pion 
scattering deserves mention. The Roy equations 82 have been reanalyzed in great detail in Ref. 
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A similar analysis was also performed in [51] . The results of this analysis have been combined with the 
results of the order p 6 calculation mentioned here to get a series of very precise predictions for the pion- 
pion scattering system in [53] EH] ■ These predictions have been nicely confirmed by experiment [571 155] 
and [89J. The relevance for the mechanism of spontaneous symmetry breaking in QCD is discussed in 
|55j . The case with a possible small value of the quark condensate jHj is now ruled out. 

The work of [BUI [5J3 [HE] has been criticized in Refs. [HHl US] on two grounds, the high-energy input 
used as well as the fact that certain sum- rules were not obeyed by the results of [52] [53] [5^]. The 
criticisms in these papers have been answered in Refs. |92[ I93|. They showed that the high-energy input 
used in [53 [5J3 EH] also satisfied the requirements as well as the one advocated in [Hill Ell- They showed 
as well that changing the high-energy input to the one of jHO] IHIj did not change their result outside 
quoted errors. The sum-rules discussed in [Hill HH] are not well convergent. Ref. [521 EH] showed that 
more convergent versions of the sum-rules are well satisfied by their results. 

The successful prediction of pion-pion scattering at low-energy from the combination of ChPT and 
Roy equations is one of the great successes of theory in low-energy hadronic physics. 

3.6 Pion form-factors 

The pion vector and scalar form-factors are also known analytically to order p 6 [73]. They are defined 
respectively by 

(n\p 2 )\uu + dd\^( Pl )) = ^F s (s n ), 
(7i i (p 2 )\-(u^u-d-f fl d)\TT j ( Pl )) = ie i3 \p 1 ^ + p 2 ^F v (s n ) , (108) 

where s n = {p2—pi) 2 - The scalar form-factor is defined with an isospin-zero scalar source. The isospin- 
one scalar form-factor can be defined analogously but it only starts at 0(p 4 ). The vector form factor 
is here defined as the isovector part only, and what we calculate here is its I z — component. Similar 
definitions exist for the other isospin components. In the limit of conserved isospin, these components 
are the relevant ones. The tree level results have been long known and it is one for the vector form- 
factor. The order p 4 result has been derived in Ref. [2] and the full order p 6 expressions have been 
obtained by the authors of Ref. [73]. Similar to the case of pion-pion scattering, the expressions are 
fully given in terms of elementary functions. 

The two form-factors have a different behaviour. The vector form-factor has relatively small correc- 
tions coming from the loop diagrams and has the corrections dominated by the part given by the LECs 
of order p 4 and p 6 . The fit to data is dominated by the spacelike measurements of the pion form- factor 
at CERN by NA7. In Fig. ^2 we have shown the quality of the fit achieved. 

This fit also allowed for a good measurement of a p 4 and p 6 LEC. The result is |ZS] 

l 6 = 16.0 ±0.5 ±0.7, 
r r V2 (M p ) = -4d 51 +4c r 53 = (1.6 ± 0.5) ■ 1(T 4 . (109) 

For the pion scalar form-factor there is no direct experimental information. The momentum depen- 
dence can be derived using dispersive methods [M]. More recent determinations and discussion using 
the same method can be found in PJ] EH EJJ . The result for the p 4 LEC that follows from a reasonable 
range for the scalar radius is 

[ 4 = 4.4 ±0.3 (110) 

if one assumes a range for the scalar radius 

(r 2 y s = 0.60 ± 0.03 fm 2 . (Ill) 
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Figure 11: The spacelike vector form-factor of the charged pion and the best fit of the ChPT 
order p 6 expression enhanced by a cubic analytic term. Figure from Ref. j73| . 



An approximate value for an order p 6 constant could also be obtained from this comparison, but its 
actual value was dependent on the other p 4 LECs input values used. Ref. |73] obtained 
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The last full calculation in two-flavour ChPT to order p 6 I am aware of is for the radiative decay of the 
pion in Ref. [HE]- The process ir + — > £ + i"y at lowest order is nothing but the QED correction to the 
pointlike ir — > iv decay. First at order p 4 are there contributions from the pion structure and the allow 
a sizable effect for tt + — > e + vy since the QED Bremsstrahlung contribution is helicity suppressed there. 
The order p 4 contribution has been calculated in Ref. [2]. The order p 6 calculation was added in 
Ref. 98J. Again, a rather good convergence from the order p 4 to the order p 6 result was seen. The 
combination of order p 4 LECs that can be obtained from this calculation is |9*8] 



2u — Ir 



0.00315 ±0.00030. 



(113) 



This was derived from the value for the axial form-factor in the decay. There have since been new 
data from the PIBETA collaboration [99 which are compatible with the value of the axial form-factor 
used in [73] to determine (|113|) . Note that the general agreement of the data with the distributions 
predicted by ChPT, or the standard V — A picture is not very good. The variation of the form-factors 
with momenta is expected to be small and the experimental variation is significant, alternatively other 
form-factors can contribute. This discrepancy is at the edge of statistical significance as discussed in 



3.8 Values of the low-energy constants 

The LECs of order p 4 were first determined in the original paper [2] using the then best values and 
the order p 4 expressions that were derived there as well. All these quantities are now known to order 
p 6 which in principle allows for a much more precise determination. A problem that surfaces at this 
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level is how to deal with the values of the unknown order p 6 LECs. In all of the papers mentioned 
earlier similar estimates using resonance saturation have been used. Since quark-mass corrections in 
the two-flavour case are suppressed by powers of M% the unknown parameters have typically a fairly 
small effect but no general study of this has been undertaken. 

The best estimate of l\ to l± comes from the combination of ChPT at order p 6 and the Roy equation 
analysis of Ref . [HE] • The other two known parameters are I5 and Iq as discussed earlier. The best values 
at present are thus 





ll = 


-0.4 ±0.6, 




h = 


4.3 ±0.1, 




h = 


2.9 ±2.4, 




h = 


4.4 ±0.2, 


fe- 


-h = 


3.0 ±0.3, 




k = 


16.0 ±0.5 ±0.7. 



(114) 

The value for T 4 is in complete agreement with the one in Eq. (|110j) . The error is smaller due to the 
fact that the work of [HE] allowed to pin down some of the input better than in [73]. The remaining 
parameter I7 is not known directly from phenomenology. It can in principle be determined from lattice 
QCD evaluations of the neutral pion mass in the presence of isospin breaking. Its value was estimated 
to be about 

Z 7 ~5-10" 3 (115) 

from 7r°-?7 mixing in j2]. 

Some combinations of order p 6 LECs are known as well. Two have been derived from the vector and 
scalar form- factor of the pion and have been given above. Ref. jHHl has derived some more combinations 
which are implicit in the values of 65 and b$ quoted there. 



4 Three-flavour calculations 

The three-flavour calculations started soon after the first two-flavour results had appeared. In the two- 
flavour case all calculations have used the same renormalization procedure as discussed in Sect. 12.101 
with the ChPT modified version of modified minimal subtraction. This is unfortunately not true for the 
three-flavour case where at least three different renormalization schemes have been used. These schemes 
are related and can in principle be related to each other. In practice it has made use of calculations 
of different groups of authors difficult. The situation at present is that most calculations performed by 
Bijnens and collaborators have been used together with a similar scheme of analysis and the standard 
ChPT subtraction scheme. Most of the other calculations have done only a rather small amount of 
numerical analysis, often using a set of LECs determined at order p 4 . 

In this section I give a list of the calculations which have been done and show a few typical numerical 
results. The section finishes with a discussion of how the order p 6 LECs have been treated in the existing 
calculations and I also briefly discuss more recent developments which have happened in this regard. 

4.1 The vector two-point functions 

The vector currents are defined as 

W-^', (116) 
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where the indices i and j run over the three light quark flavours, u, d and s. Working in the isospin 
limit all SU(3) currents can be constructed using isospin relations from 

v;{x) = l=(y}\x)-V™(x: 



TO = -4W 1 (x) + Vf(x)-2Vj-i.r 



r33< 

V6 V> V ^ 7 ' '" V ^ 7 ~ V 



vf(x) = vf(*). (ii7) 

These will be referred to as the isospin or pion, hypercharge or eta and kaon vector currents respectively. 
All others can be defined from these in the isospin limit. E.g. the electromagnetic current corresponds 
to 

vr = -jf» + -jf» • ( iig ) 

The two-point functions are defined in terms of the currents as 

n£(g) = i\d'x e** (0\T(V^x)V u a (0)y\0) , (119) 

for a = ir, f], K . All other vector two-point functions can be constructed from these using isospin 
relations. Lorentz-invariance allows to express them in a transverse, IT^, and a longitudinal, n^, part 
via 

K a = (» - sV)n&(g 2 ) + q,q^vl(q 2 ) • (120) 

For a conserved current, which is the case for the pion and eta current, the longitudinal part vanishes 
because of the Ward identities. 

These two-point functions were first evaluated in Ref. |100| for the isospin and hypercharge, i.e. the 
diagonal, vector-currents. The missing case for the vector-current with Kaon quantum numbers was 
later added by Ref. |l()lj and |1()2| . The relevant diagrams to order p 6 are shown in Fig. E3 There is 
no order p 2 contribution, the order p 4 diagrams are the first line, (a-c), and the remainder are the order 
p 6 diagrams. 

There are no proper two-loop integrals needed to evaluate the vector two-point functions to order 
p 6 . They have not been used very much for phenomenological purposes but could in principle be used 
as low-energy constraints on sum-rule analyses in this channel. Some results were given in Ref. [100J. 
The papers |1U1J and [100. I1U2"] used a different subtraction scheme so a full comparison has not been 
done, but the checked parts are in agreement between the two independent evaluations. 

The isospin breaking vector two-point function was calculated to order p 6 in Ref. |lU3j and were 
used as constraints in a sum- rule analysis of isospin breaking in vector meson decay constants J1U4J . 

4.2 Scalar two-point functions 

The scalar two-point functions can be calculated from a similar set of diagrams as the vector two-point 
functions by replacing the insertion of the vector currents in Fig. ^] by scalar currents. The main one 
used in this sector is 

U s (q)=i j d A xe iq - x (0\T((S uu (x) + S dd (x))S ss (0y\0) , (121) 

with the scalar densities S 1 ^ defined in Eq. (|12|). It was calculated to order p & by Moussallam |lU5j and 
used in an analysis to obtain bounds on L\. 
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Figure 12: The Feynman diagrams needed to calculate the vector two-point functions at 
order p 6 . The crossed circle stands for the order p A vertex insertion. Wiggly lines are the 
external vector currents. Dots are order p 2 vertices and a square is an order p 6 vertex. The 
solid lines are meson propagators. Figure from Ref. |1U1J . 



4.3 Quark condensates 

Another quantity which can be evaluated without proper two-loop integrals at order p 6 is the quark 
condensate 

(*V>. (122) 

It has been evaluated in the isospin limit in Ref. [106J and the isospin breaking corrections in Ref. |lU7j 
for the three possible cases i — j — u,d, s. 

The value of the quark condensate depends on the constant H[ at order p 4 . This cannot be directly 
measured in any physical process. It's value depends on the precise definition of the quark densities in 
QCD. To order p 6 the local counterterms of that order also contribute to (c^qi). The corrections when 
L\ and L r & are assumed to be zero are fairly small as can be seen from the plots in Ref. [106J for (uu) 
and (dd) but they were sizable for (ss). 

As discussed shortly in Sect. El it is now clear that the quark condensate in the two-flavour case 
remains large also in the chiral limit where m u and rrid are sent to zero. The equivalent question 
for the three-flavour chiral limit remains open as discussed in jHS| and [HH1 1108] • The situation here 
is analogous to the situation in the two-flavour case before the latest results on pion-pion scattering 
threshold parameters, all results seem to indicate that the standard picture as described in this review 
is consistent but an alternative scenario is not ruled out. The main remaining obstacle in the three 
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flavour case is that the effects of the order p 6 LECs have not been constrained in sufficient detail yet to 
see how this question can be resolved. Work is ongoing whether the results for ttti and irK scattering 
can be sufficiently refined theoretically to provide such a test. 

4.4 Axial-vector two-point functions, masses and decay- constants 

These are the simplest quantities requiring proper two-loop integrals. The axial-vector currents are 
defined as 

A%(x)=tf 7(t75 q1 (123) 

where the indices i and j run over the three light quark flavours, u, d and s. Working in the isospin 
limit all SU(3) currents can be constructed using isospin relations from 



a;(x) - 


= ^(W-WJ 
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Al(x) - 


= ^{Vi\x) + V?{x)- 


-2V™(x 


<(*) = 


- V^(x). 





(124) 

These will be referred to as the isospin or pion, hypercharge or eta and kaon axial-vector currents 
respectively. All others can be defined from these in the isospin limit. 
The two-point functions are defined in terms of the currents as 

KM = *Jd A x e*** (0|T(^(x)^(0))t|0> , (125) 

for a = 7T, rj, K. Using isospin relations, all other axial- vector two-point functions can be constructed 
from these using isospin relations. Lorentz-invariance allows to express them in a transverse, iT 1 ), and 
a longitudinal, IT^, part 

K = («Wfc - <? V)n2(g 2 ) + fcfcJiiEfo 8 ) • (126) 

The axial currents are only conserved in the chiral limit, i.e. when the masses of all quarks in the 
currents involved vanish. Thus in general there is both a longitudinal and a transverse part. 

The axial-vector two-point functions have contributions from one-particle-reducible diagrams, those 
where the cutting of one line makes the diagram become disconnected. The full set is shown in Fig. ED 

The filled circles in Fig. [EI] are the sum of all 1PI diagrams. Those up to order p 6 are shown in 
Fig. EJ The diagram in Fig. ITlT k). the so-called sunset diagram is the first diagram we encounter 
in this section that requires the evaluation of proper two-loop integrals. Methods to perform these 
integrals thus needed to be developed. The methods derived in [HE] need to be generalized to the case 
with different masses in the loop. An efficient method to get the sunset diagram at zero momentum 
was derived in Ref. |lU9j using recursion relations between the various integrals. A different derivation 
of the same relations was presented in Ref. [101J. The method for the momentum dependence of the 
sunsetintegrals of jHS] was extended to the different mass case in |ll)lj . Ref. |ll(Jj used a different 
variation on the sunset integral method of JSH] to perform their numerical analysis. 

In Ref. [110J the axial- vector two-point function with pion and eta quantum numbers was evaluated 
to order p 6 . They also used this to evaluate the pion and eta masses and decay constants as well as a 
simple sum-rule analysis [111 . The axial- vector two-point functions were calculated also in Ref. |1U1J 
where the same quantities with kaonic quantum numbers were also evaluated to order p 6 . At this level 
all masses and decay constants were known in the isospin limit. Note that the masses can be calculated 
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Figure 13: The diagrams contributing to the axial-vector two-point function. The filled 
circles indicate the one-particle- irreducible (1PI) diagrams. Solid lines are pseudoscalar 
meson propagators and the wiggly lines indicate insertions of an axial- vector current. For 
the inverse propagator the wiggly lines are meson legs and for the decay constant the right 
wiggly line is a meson leg while the left remains an axial current. Figure from Ref. |1U1J . 
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Figure 14: The set of diagrams contributing to the 1PI quantities, (a) axial-vector-axial- 
vector (b) axial- vector-pseudoscalar (c) pseudoscalar-pseudoscalar. (d)-(l) the respective 
diagrams when the dashed lines are replaced with the external legs of (a), (b) or (c). A line 
is a meson propagator, a wiggly line an external source, a dot a vertex of order p 2 , a crossed 
circle a vertex of order p 4 and a crossed box a vertex of order p 6 . Figure from Ref. [101] . 
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from the position of the pole in the full propagator and the decay constant by direct evaluation of the 
matrix element 

{Q\Al\M{p))=iy/2F M P lt . (127) 

Flj can be determined as well from the residue of the meson pole in the longitudinal part of the axial- 
vector two-point function. That all these methods of calculating the masses and the decay constants 
give the same answer was explicitly checked in Ref. |1U1J . 

What was found for the masses was that the order p 4 corrections were small for the standard set 
of input parameters from fits at order p 4 , see Ref. |112j for this determination. However, the order 
p 6 corrections were rather large for both |ll()j and |l()fij . It should be noted that there was analytical 
agreement between those two references for the parts that could be checked without relating the different 
renormalization schemes and the different way of evaluating the integrals. At present there is no obvious 
solution for the presence of these large corrections. It is possible to choose the order p e LECs to cancel 
these corrections but the LECs involved are those coming from the scalar sector and are the most 
difficult ones to estimate. The naive estimates used in |l(Jlj gave a very large range for the contribution 
from the order p 6 LECs, from extremely large to zero. In Sect. 14.81 1 will discuss a few more relevant 
results. 

The masses and decay constants away from the isospin limit have also been worked out. These 
calculations were reported in Ref. |!U7j . 

4.5 K M 

The decay K — ► Tnriu was first treated using current algebra methods to obtain the lowest order 
result. The order p 4 calculation was performed by two groups simultaneously |113[I114] with complete 
analytical agreement. The reason for calculating this decay beyond lowest order was two-fold. The 
most accurate determinations of the form-factors showed a significant deviation from the lowest order 
prediction. It was expected that the order p 4 calculation would allow to fit the experimental results 
and allow for a more accurate determination of L\ with % = 1,2,3. 

The corrections to the lowest order were of the expected order but not very small. This prompted an 
investigation of higher orders using dispersive methods. This was done in Ref. |115j . In that reference 
also the missing form-factor R defined below was evaluated to order p 4 . The higher order corrections 
were of the expected order indicating a converging series for this decay. The fit results of Ref. |115j 
have been the standard set of values for the order p 4 LECs replacing the earlier full fit of Ref. [3]. 

When the full work on pion-pion scattering in two-flavour ChPT was finished it became interesting 
to check whether this calculation was compatible with the values of the order p 4 LECs determined from 
K14. To be able to do this K ti needed to be determined also to order p 6 . This calculation was performed 
by the authors of Ref. ilOol lllb] . I will now present the results from these references. 

The K14 decay processes are 

#+(p) _> n+(p + )n-(p_)£ + (p e )v e (p v ), 

k + (jp) - K°( P+ )K (p-)£ + (PiMp»), 

K\p) -► n-(p + )n°(p^£ + (p e Mp u ). (128) 

The corresponding momenta are given inside the brackets after each particle. Following the original 

work |117j . the K^ decays are parameterized in terms of five kinematical variables: 

i) Sn, the squared effective mass of the dipion system. 

ii) S£, the squared effective mass of the dilepton system. 

iii) 8^, the angle between the ir + and the dipion line of flight with respect to the Kaon rest frame. 

iv) 0£, the angle between the £ + and the dilepton line of flight with respect to the Kaon rest frame. 

v) 0, the angle between the n — tc and e — v planes with respect to the Kaon rest frame. 
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However, from the point of view of the hadronic system alone, a better choice of variables is to use 
in addition to s ff , cos 6^ and (ft, also 



U = (p+ ~ pf and u^ = (p_ - p) 2 , (129) 



related through 
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(130) 



O 7T \/ ■'- l ' ' ' -.- / ^TT 5 



y/l - 4m2/s Tl 
A(mi,m 2 ,m 3 ) = m^ + m^ + rn\ — 2 (r7i 1 m 2 + m 1 m 3 + m 2 m 3 ) . (131) 



With the previous notation the amplitude for the decay K + — > 7r + 7r £ + z/<? is 



with 



^V: s u(p v ) ltl (l - l5 )v(Pi)(V" - A") , (132) 



V^ = 3-e M! ,pa(P£ + 2V)"(P+ + P-) P b+ - P-Y , 

^ = -— [(p + +p-V J F+(p + -p-) M G + (p / +p„)^ J R]. (133) 

% 

Vu S is the relevant CKM matrix-element. The other two amplitudes, T~° and T 00 , are defined similarly 
Here we are interested in the F and G form-factors. The H form- factor is known to 0(p 6 ) |118j and 

the R form-factor always appears with a factor of m\ and its contribution negligible for the electron case. 

R is known to order p A |115j . The form-factors are functions of s n , Si and cos9 n only, or alternatively 

of s nj tjr and u^. 

The relations between the form- factors and the intensities, easier to obtain from the experiment, 

can be found in |119j or in |115j . 

The amplitudes for the three processes of Eq. (J128J1 are related using isospin by 

T +- = ^_ + T oo ( 134 ) 

v2 

with T* J the matrix element defined in Eq. (J132J) . T~° is anti-symmetric under the interchange of the 
pion momenta while T 00 is symmetric. This also implies relations between the form-factors themselves. 
Observe the different phase convention in the isospin states compared to |115j where M 00 appears with 
a minus sign because of the Condon-Shortley phase convention. 

The form-factors F and G can be decomposed in partial waves. The partial wave expansion is not 
simply for F and G since the components with a well defined L z = 0, ±1 need to be expanded |117[I119| . 
The relevant expressions can be found in |117t I119| I115| I10(j| and a discussion of their relative sizes and 
how they can be used in experiment is in the mentioned references but more elaborated in Ref. 
I will not go further into this but simply quote the experimental results. 
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Figure 15: (a) One-particle irreducible tree level diagram, (b) One-particle irreducible one- 
loop diagrams. Dots refer to strong vertices or current insertions from £ 2 > £4 or £ 6 . External 
legs stand for pseudoscalar or weak current. Internal lines are pseudoscalars only. Figure 
from Ref. fTUfij . 






<g) 





(hi 



Figure 16: One-particle irreducible order p 6 diagrams with only one- loop integrals. The 
dots are order p 2 vertices except the top dot in (d) and (e) which is an order p 4 vertex. In 
addition there are the diagrams (b) and (c) of Fig. [T3] with one of the dots replaced by a 
0{p i ) vertex. Figure from Ref. [106J. 



The lowest order calculation was done using current algebra methods by Weinberg |121j . Both the 
F and G form-factor are equal and given by a single insertion of a £2 vertex in diagram (a) of Fig. [TBI 
The result is 

F = G = 4=^, (135) 

where tuk is the physical Kaon mass. At higher orders more diagrams become relevant. The tree- level 
and one-loop diagrams are shown in Fig. E3 To order p A we need in addition the diagram of Fig. flBT a) 
but with an order p A vertex as well as the loop diagrams (b,c) of the same figure with only order p 2 
vertices. With order p 4 result, the experimental data at he time can be nicely fitted, as discussed 
extensively in Ref. |114j and later in Ref. |115j . 

The diagrams at order p 6 that are relevant for the calculation of the F and G form-factors defined 
in Eq. 11331 are shown in Figs. EI] and EH in addition to diagrams in Fig. E3 with insertions of an order 
p 6 vertex in diagram (a) or one order p A vertex in (b,c). The diagrams shown in Fig. E3 do not require 
any proper two-loop integrals. They can be calculated using straightforward methods and only present 
difficulties due to the length of the expressions. 

Proper two-loop integrals show up in the diagrams shown in Fig. El The left one, (a), is called 
the sunset- diagram and the integrals needed are the sunsetintegrals. These were discussed briefly in 
Sect. 14.41 The diagram on the right-hand-side, (b), in Fig. El is the single most difficult diagram in 
this calculation. It is called the vertex-diagram or sometimes the fish-diagram. There have been quite 
a few different methods used to evaluate this type in integrals in order p 6 calculations. The one used 
in most calculations is based on the work of Ghinculov, van der Bij and Yao |122j and can be found in 
those references or in Ref. |123| . 

In addition to the diagrams shown in Figs. IToTfrTl there is the contribution from wave-function renor- 
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Figure 17: One-particle irreducible order p 6 diagrams with irreducible two- loop integrals. 
The dots are order p 2 vertices. Figure from Ref. 106 . 
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Figure 18: Comparison of the Omnes improved estimate of Ref. |115j with the full order 
p 6 of calculation of Ref. [106^ using the same set of values of L\ as input together with 
Fff = 93.2 MeV. The shaded band is the experimental result of [SI]. Figure from Ref. 106 . 



malization from the diagrams shown in Fig. EJ The order p 6 expressions are very long and are partly 
reported in a numerical parameterization only in Ref. [106J. It was found that the parameterization 
based on a truncated partial wave expansion as proposed in Ref. |120j worked very well. 

To show the results, we first compare the Omnes improved order p 4 result of Ref. |115j together with 
the data of Ref. jHTj in Fig. E3 The result was fitted to that data so there is quite good agreement. 
However, when we compare the full order p 6 calculation we notice that there is quite a difference with 
the dispersive, i.e. the Omnes improved, estimate. We also see a reasonable convergence of F going 
from order p 2 to p A to p 6 . 

After performing a fit to the data of Rosselet et al. [HI], the order p 6 result fits the data nicely. 
The resulting different orders are shown in Fig. EH For both form-factors there is a good convergence. 
The estimate of the contributions of the order p 6 Lagrangian, the curve labeled C^-only, is also shown. 
It is small for both form-factors. The estimate for G is somewhat larger, but follows from measured 
resonance decays so it is fairly certain. The full description of how this estimate has been done can be 
found in Refs. [TmiP0T)j . 

Another feature visible in the figures is that for F, which is dominated by the S"-wave, final state 
rescattering is rather important since the contribution from the imaginary part to \F\ is quite visible. 
In contrast, G is dominated by the P-wave and final state rescattering effects are fairly small. 

After the work of Ref. |106j was finished, new data have appeared. The main new data from the 
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Figure 19: The F (a) and G (b) form-factor at sg = and cos^ = 0. The shaded band 
is the result of 81J. Shown are the full result for the absolute value, and the real part to 
lowest-order, 0(p 4 ) and 0(p e ). We also show the contribution from the 0(p 6 ) Lagrangian; 
this is dominated by the vector contribution. The curves for Re G and \G\ are very close 
since Im G is very small. Figure from Ref. [106 . 



BNL experiment [S3 EH] were included in the full fit done in Ref. [1.07J. These data were preliminary 
at the time, so they were not called the main fit in that reference. The best fit, including the newer 
data, is called "fit 10" in Ref. |lU7j and is the one which has been used in most works following it. 

In Ref. [106J, a study of the dependence on the constants L\ and L r & was done as well. A main goal 
was to see whether good fits could be obtained with these at values different from zero as well. It was 
found that such is indeed the case, and that in order to obtain these two order p A constants, additional 
experimental input was needed. This will be discussed in the context of the scalar form-factors, pion- 
pion and pion-kaon scattering below. 

We also expect that new data on various K^ decays will become available from the NA48 experiment. 



4.6 Vector or electromagnetic form-factors 

The most general structure for the on-shell pseudoscalar-pseudoscalar-vector Green function is dic- 
tated by Lorentz invariance. With the additional use of charge conjugation and electromagnetic gauge 
invariance one can parameterize the pion and kaon electromagnetic matrix elements as 

(7r + (g)|j»)> = (q»+p»)F*(t), 

(K+(q)\ Jtl \K + (p)) = (q»+p»)F« + (t), 

(K%q)\j,\K (p)) = (q,+p,)Ff(t), 



(136) 



with t = (q — p) 2 . The current j\i refers to the electromagnetic current of the light flavours 
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loops only L[ = C[ = 
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Figure 20: The real and imaginary parts of the loop diagrams at order p e and order p 6 with 
all L\ — and CI = 0, for the pion form-factor. Notice that there is convergence both for 
the real and the imaginary part. Figure from Ref. J123J . 



The quantities Fy, Fy and Fy will be referred to hereafter as the vector form- factors or simply the 
form-factors. They are also defined in the crossed channel {0\j fl \M a (p)M b (— q)). 

To lowest order the vector form-factors are constant and are simply the charge of the relevant 
meson under the current under consideration. The order p A calculation was performed by Gasser and 
Leutwyler in Ref. |124j and confirmed for the pion form- factor in Ref. [HS1 • The full order p 6 calculation 
was performed independently by Bijnens and Talavera |123| and Post and Schilcher |125| 11261 D-27J. 
The results of Post and Schilcher for the neutral kaon form-factor are in Ref. |126j and for the pion 
form-factor in the appendix of Ref. [127J. The radius was also derived for the charge radius of the 
combination of vector form-factors which has only higher order breaking in the quark-masses as derived 
by Sirlin |128[I129T| in Ref. [125J. The two calculations use a very different method to calculate the sunset 
and vertex integrals. In addition, Post and Schilcher have used a slightly different subtraction scheme 
as well as only older values of the order p 4 LECs. The analytical results which could be compared 
between the two calculations agree. 

The numerical results for the pion vector form-factor are very similar to those of the two-flavour 
calculation discussed in Sect. 13.61 The contribution from the pure loop diagrams of order p 4 and order 
p 6 is shown in Fig. 1201 They are rather small at each order and the convergence is nice. The order p 6 
contribution is significantly smaller than the order p 4 result. 

The calculation can now be used to fit to the available data. The fit is very similar to the one in the 
two-flavour case shown in Fig. ^2 The data included in the best fit as shown in Fig. |^ The inclusion 
of time-like data does not really change the fits as discussed in [73] and |123j . Figures for that part of 
the fits can be found in those references. Note there seems to be some numerical problems in the work 
of Ref. |127j . The real part of the result is compared with the measurements in the appendix of that 
reference but the imaginary part quoted there is much larger than thee one found in |123j and way too 
large to be compatible with the known pion phase-shift at that order in ChPT. 



The fits of Ref. 



allow to determine the parameter L r 9 to order p as 
L£(0.77 GeV) = (5.93 ± 0.43) x 10~ 3 . 



:i38) 



The error due to experiment only is about half this, the remainder is mainly from the estimate of the 
order p 6 parameters that contributes proportional to m^ or m 2 K . The curvature also has a contribution 
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Figure 21: The comparison of the space-like measurements with the ChPT calculation for 
the pion electromagnetic form-factor. Notice that there is excellent convergence over the 
whole kinematical range. The curve labeled p 8 includes a term Cft 3 in the fit in addition to 
the full order p 6 expression. The data are from Refs. J13U|. I13l| . Figure from Ref. 



from an order p 6 constant allowing a direct determination giving [123]. 
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(0.22 ±0.02) x 10- 



(139) 



This result is in good agreement with the various resonance estimates discussed below. Eqs. (J138J) and 
(I139J) are the direct three-flavour equivalents of the two-flavour results of Ref. [ZH] given in Eq. ()109|) . 

The kaon electromagnetic form-factors have also been calculated. The data here are not so good and 
there seem to be some problems with the absolute normalization. The data from Refs. [1311 1132[ I133J 
are in agreement with the results from ChPT within the experimental errors. 

Electromagnetic form-factors are near t = often described by the charge radius defined by 



(r 2 Yv K+ ' K ° 



b dt V 



(140) 



t=o 



and the data and fits given can be used to extract /predict the various radii. This has been done in 
Refs. jZHl H23j for the two- and three-flavour case respectively. 

The K° electromagnetic radius has been measured in a kaon regeneration on electrons experi- 
ment |134j with the result 

(r 2 )f = (-0.054 ± 0.026) fm 2 . (141) 

The decay K^ — > n + TT~e + e~ also has contributions that contain the neutral kaon electromagnetic form- 
factor but the extraction from the data is rather model dependent and I will not use these data. This 
radius was predicted in Refs. |126| I123J , the result from |126j used older values of the order p 4 LECs. 
The prediction from |123j is 
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(r 2 )f 



{ - 0.0365 [loops p A ] - 0.0057[loops p 6 ]} fm 2 
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Ft 



r>K° 
rL V2 



-0.042 ± 0.012) fm 2 



(142) 



The error is based on assuming that the unknown contribution is not larger than twice the order p 6 
loop contribution. The result is in good agreement with the measurement ()141|) . We can also turn the 
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argument around and obtain 

Rf 2 = (-0.4 ± 0.9) x 10" 5 GeV 2 . (143) 

This result is in reasonable agreement with the resonance estimate presented in Ref. |123j . 

4.7 K G 

The decays considered in this subsection are 

tf+(p) _> 7rVK + (wH(Pv) , [K& (144) 

lfO(p) _ 7r-(p0^ + ^)^(p,) , [Kl\ (145) 

and their charge conjugate modes. £ stands for \x or e. The short notation for each decay is given in 
the square brackets. 

The matrix-element for K^, neglecting scalar and tensor contributions, has the structure 

T = ^KWV,p)> (146) 

with 

F = u(p v )^(l- l5 )v( Pe ), 
F, + ( P ',P) = <n°(p')\V^ 5 (0)\K + (p)>, 

= -^i(p+p)»ff 7T °(t) + (p-p') fl f K+ « (t)}. (147) 

To obtain the K® 3 matrix-element, one replaces F+ by 

F,°(p',p) = <n-(p')\V^(0)\K°(p)> 

= (p'+p)J? 7r ~(t) + (p-p>) IM f I f--(t). (148) 

The processes ()144|) and ()145|) thus involve the four Kg 3 form-factors f± n (£), f± n (t) which depend 

on 



t = (p'-pf = (Pz + Pv)\ (149) 



the square of the four momentum transfer to the leptons. 

Only the isospin conserving part is known to order p 6 . The isospin breaking corrections are one of 
the remaining unknowns in the determination of the CKM matrix-element V us . In the isospin limit we 
have that 

f± = ft = f± +7T ° = /±° *~ ■ (150) 

jKir j g re f errec [ t as t^ vector form-factor, because it specifies the P-wave projection of the crossed 
channel matrix-elements < | V^~ l5 (0) \ K + , 7r° in >. The S"-wave projection is described by the scalar 
form-factor 

/oft) =/+(*)+ 2 t 2 /-(*). (151) 

Analyses of K^ data frequently assume a linear dependence 



f + , Q (t) = / + (0) 



1 + ^+,0 o~ 

mt. 



(152) 
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For an early discussion of the validity of this approximation see |3] and references cited therein. As 
pointed out in Ref. \M$ and checked afterwards by the ISTRA+ [U)], KTeV HHZj and NA48 [T3K] . 
the linear approximation is not sufficient at the present level of precision. 

The form-factors f±,o{t) are analytic functions in the complex t-plane cut along the positive real 
axis. The cut starts at t = (mx + rn n ) 2 . In the phase convention used here, the form-factors are real in 
the physical region 

mj < t < {vfiK — m n ) 2 . (153) 

A discussion of the kinematics in K^ decays can be found in |112j and references cited therein. 
The total result can be split by chiral order 

m = f?\t) + f?\t) + f?\t)i (< = +,-, 0). (154) 

The lowest order result has been known for a very long time and is fully determined by gauge 
invariance. 

/| 2) (t) = / (2) (t) = l,, /i 2) (t) = . (155) 

The order p A contribution was first calculated within the ChPT framework by Gasser and Leutwyler |124| 
with earlier results by Leutwyler and Roos |139| . The result contains the nonanalytic dependence in the 
symmetry parameters predicted by |14Uj . Partial studies at order p 6 have also been done, the double 
logarithm contribution is small as was shown in Ref. [HE] and a possibly large role for terms with two 
powers of quark masses has been argued for in Ref. |141j . The latter reference also pointed out the 
strong interdependence of the measurements of Fk and V us as a possible solution to the CKM unitarity 
problem. 

There exist two full calculations at order p 6 . The work of Post and Schilcher |127j and Bijnens and 
Talavera |135| . The first work uses outdated values of the ChPT constants as well as an older version of 
the classification of p 6 constants. Due to the different methods of subtraction and splitting up the various 
proper two-loop integrals in an analytical and numerical part, a full analytical comparison between the 
two results has not been done. The parts which can be compared are in agreement analytically. The 
numerical results of both papers are in agreement for / + when the same input values for all parameters 
are used. There is a disagreement for /_ but as discussed in Sect. 14.61 there are some indications for 
numerical errors in Ref. |127j . 

Let me now discuss the main results of Ref. |135j . I will present some of their numerical results later 
but first discuss the mo del- independent relations at order p G . There are a large number of order p 6 LECs 
contributing to both form-factors, their contributions to the processes discussed in the previous, this one 
and the next subsection satisfy several relations. The combination which is relevant for the curvature 
of f+(t) is the same combination that appears in the curvature of the pion electromagnetic form-factor. 
Its value has been determined from the data and is given in Eq. ([139)1 . That way Ref. |135| predicted 
the curvature in f+(t). In a similar way, the curvature of the scalar form-factor, fo(t), is determined by 
the order p 6 constants C[ 2 . Its value can in principle also be determined from the curvature of the pion 
scalar form-factor [142J, Sect. 14.81 There are also relations between the order p 6 constants contributing 
to the various charge radii. This relation is essentially the Sirlin relation of |128M129] . 

A more surprising relation that connects the values of /+(0) with the slope and curvature of the 
scalar form- factor fo(i), was discovered in Ref. |135j . One constructs the quantity 

/o(*) = /+(«) + 2 t 2 (/-(*) + 1 - Fk/F*) = f Q (t) + 2 t 2 (1 - F K /F V ) . (156) 

This has no dependence on the L[ at order p 4 , only via order p 6 contributions. Inspection of the 
dependence on the C\ shows that 

fo(t) = 1 -^j(C[ 2 + C^(m^- mlf + 8±(2Cl 2 + C r u )(m 2 K + mj) 

■K TV 
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— t 2 C[ 2 + A(t) + A(0). (157) 

TV 



It should be emphasized that the quantities A(t) and A(0) can in principle be calculated to order p 6 
accuracy with knowledge of the L\ to order p 4 accuracy. In practice, since a p 4 fit will include in the 
values of the L[ effects that come from the p 6 loops (due to the fitting to experimental values) we 
consider the p 6 fits to be the relevant ones to avoid double counting effects. 

The definition in (J156J1 has essentially used the Dashen-Weinstein relation |143j to remove the L\ 
dependence at order p 4 . It has also the effect that it removed many of the C\ from the scalar form- factor 
as well. The corrections which appear in the Dashen-Weinstein relation are include in the functions 
A(t) and A(0), these have both order p 4 |124t 1140] and order p 6 contributions. 

It is obvious from Eq. (J157J) that the needed combination of CI can be determined from the slope 
and the curvature of the scalar form-factor in K^ decays. 

It seems possible that C[ 2 can be measured from the curvature of the pion scalar form-factor near 
|142j . With this calculation is complete, one can use the dispersive estimates of the pion scalar form- 
factor together with only a Ao measurement in K^ to obtain the p 6 value for /+(0). There are also 
some dispersive estimates for the relevant scalar form-factor. Unfortunately, these were not in a usable 
form |144j when Ref . |135j appeared, but have since been treated in Ref . |145| . Other estimates of the 
relevant constants are the one used in |135j coming from Ref. |139j and the resonance chiral theory 
result of 1461. 



One feature that is visible in Eq. 11571 is that the value of /+(0) only differs from 1 by terms of order 
(m 2 K — m 2 ) 2 . This is not only true for the analytic contributions written out explicitly but also for the 
entire expressions for /+(0). This is known as the Ademollo-Gatto theorem |147j and holds to all orders 
in ChPT. 



Let us now discuss the adequacy of the linear parameterization for the form-factors. In Ref. 
it was pointed out the measured value of a form-factor at zero and the slope are rather dependent on 
the curvature, even if the data themselves do not show the presence of curvature. The fitted value of 
the slope and the form-factor at zero can change significantly outside the quoted errors due to this. 
This was shown in Ref. J135J for the case of the then best available data in both the neutral decay 
channel, CPLEAR jHH] and charged decay channel, KEK-PS E246 [H]]. This is shown in Fig. E2 for 
the KEK-PS E246 data. The discussion for the CPLEAR data can be found in Ref. |135j . The size 
of this effect is fairly small for the KEK-PS data but much larger for the CPLEAR data. It is very 
important to get data at small t to minimize this effect as far as possible. 

Since Ref. |135| appeared there have been newer high precisions data from the ISTRA+ 136J, 
KTeV |137j and NA48 |138j collaborations. The last two have clearly demonstrated the existence of 
the curvature experimentally but there are some disagreements in their results. The question above is 
relevant now at a higher level of precision. How much does the presence of a t 3 term affect the latest 
high precision fits? 



A recent reference discussing the various results and impact on the determination of V us is 
The final theoretical error has been significantly improved since the work of Leutwyler and Roos 
and is now roughly one percent. The bulk of this comes from the uncertainty in the needed order 
p 6 constants. More theoretical work on dispersion relations in the scalar form-factor together with 
experimental studies of the scalar form-factor together with the relation discussed above should allow 
to pin down this contribution to better precision. The uncertainty on the theory side is definitely smaller 
than the violation of CKM unitarity when taking V u d from neutron beta decay. 
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Figure 22: The KEK-PS E246 data together with the ChPT result, with the linear coef- 
ficient fitted and curvature predicted from ChPT as well as the linear fit of the KEK-PS 
collaboration. 



4.8 Scalar form- factors 

The scalar form-factors for the pions and kaons are defined as follows: 

(M 2 { P )\q l q \M l {q)) = F^ h {t) 

with t = (p — q) 2 and i,j = u, d, s. M 1; M 2 are meson states with the indicated momentum. 
In the case of isospin symmetry m u = rrid = in, the various pion scalar form-factors obey 



(158) 



*s® 



lKu{t) = 2*2?* (t) = 2i^ + (i) = 2Ff d *\t) 
-2v / 2F^ + (t) = 2v / 2F^"(t), 



Ss 



Tp'K TV 77I7T 1T 

SS ~ SS 



(159) 



The kaon currents are related by the rotations in flavour space 



F s K u (t) 
F s K d (t) 
(t) 

(t) 



Sd 



F K+K+(t) 
uu 

dd 



F,y*' (t) 



Fzr^it) 



dd 

F K°K° 

uu 



W K 



F«^(t) = V2F. 
FsKt) + F s K d (t) . 



(t) 
(*) 



su 



(t) 



(160) 



The other scalar form-factors can be obtained from the above using charge conjugation and time reversal. 
Eqs. (I159J1 and ()160|) also show the notation used for the form-factors in the remainder of this subsection. 
The scalar form-factor Flf n (t) is proportional to the form-factor f (t) used in K^ decays, Sect. 14. 7| 
The scalar form-factors obey a relation similar to the Sirlin [.1281 H29j relation for the vector form- 
factor 11421 



m s — m) 2 ) . (161) 



Ki 



1K1T1 



Fs(t) ~ 2F§,{t) - 2Ff d (t) + 2F&(t) - AF^{t) = O 
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The proof of this relation is in App. A of Ref. 142 and uses a method similar to the one given in 
Ref. (129j . The data at present do not allow to test this result. 

The values at zero momentum transfer are related to the derivatives of the masses w.r.t. to quark 
masses because of the Feynman-Hellman theorem (see [124J) 



F S(°) = lL<> ^ s (0) = ^-m 2 



dm nl bsK ' dm s 



w 1 



dm u K bsK ' dm s 



F su(0) = jr-ml , F*(0) = —m\ , 



F ^ = tk m «- (162) 

The masses are known to order p 6 as discussed in Sect. 14.41 and there were sizable corrections at order 
p 6 . This also shows that large corrections at t — can be expected. The argument goes as follows. If 
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B m s + (3{B m s f + -f(B m s f , (163) 



then 

Fg(0) ^B + 2/3B m s + 3 7 (A)m s ) 2 . (164) 

So we see that in the scalar form-factors the relative p 6 corrections can get enhanced by factors of order 
3 compared to the masses. 

The lowest order results have been long known. Order p 4 results were obtained in |124[ I150J . The 
actual calculations of the scalar form-factor to order p 6 in three-flavour ChPT were performed by the 
authors of Ref. |142j . As expected from the argument just presented, large corrections were found to 
several of the form-factors involved, especially for the kaon ones. The results are also quite sensitive 
to the input values used for the LECs. In Fig. I2B1 the contributions to the two pion scalar form-factors 
at order p A and p 6 are shown. The are shown for the values of fit 10, one with a changed value for 
L\ and the pure loop parts (all L[ = 0). It is clear that while for these values the corrections are not 
enormous, they nonetheless are large and show no obvious convergence in that for many t the order p 6 
contribution is larger than order p 4 . 

The analysis for the pion and kaon form-factors followed the method introduced by Ref. [01] as 
updated in Refs. |1()5| ITTHj . The actual discussion of the results can be found in Ref. [142J. I will 
restrict myself to a few small comments here. One uses the ChPT input for the form-factors at zero 
and then calculates using the dispersive methods of Muskhelishvili-Omnes its momentum dependence. 
In order to have a consistent set of other LECs when L\ and U % were varied, Ref. |142j redid the fits 
with the same assumptions as fit 10 in Ref. |107j but with range of inputs for those two L\. The form- 
factors at zero momentum were then calculated from the order p 6 ChPT expression and the momentum 
dependence calculated from those with the dispersive method using [1051102] result. The comparison of 
the dispersively estimated momentum dependence with the ChPT calculated momentum dependence 
leads to the constraint |142j 

Ll& L r 4 - 0.00035. (165) 

The comparison for the scalar radius is shown in Fig. |^] The constraint roughly describes the line 
where the ChPT and dispersive radius are in reasonable agreement. 

That the actual values of L\ and L r & has a large impact on the masses and scalar sector is shown 
in Tab. H3 The masses and decay constants for the sets of input parameters given can also be found in 
Table El The various orders in the expansion quoted there are defined as: 

(4) , e,(6) 



iWF = 1 + Fg + F 
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Figure 23: The effects of the L\ on the various contributions to (a) F<f(t)/(2Bo) and (b) 
Fg s (t)/Bo as a function of t for the case with C[ = 0. The curves are for the real parts at 
order p A and order p 6 . The curves labeled respectively Fit 10, L\ ^ and L\ = are for the 
standard values of the L\ of fit 10 in |107j . the same values but L\ = —0.003 and with all 
L\ = 0. Figure from Ref. 
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Table 3: The results for various quantities for the L[ of fit 10 and three other representative 
sets A, B and C. Adapted from Ref. [T42] . 



Set 


Fit 10 


A 


B 


C 


10 3 • L\ 




0.0 


0.4 


0.5 


0.5 


10 3 • L r 6 




0.0 


0.1 


0.1 


0.2 


FS(Q)/B (ChPT, C\ = 


0) 


2.54 


1.99 


1.75 


2.12 


FSM/Bo (ChPT, CI = 


0) 


0.020 


0.004 


-0.002 


0.008 


F«(0)/B (ChPT, CI = 
f£(0)/B (ChPT, ci = 


0) 


1.94 


1.36 


1.12 


1.51 


0) 


1.77 


1.35 


1.17 


1.45 


< r 2 >l 4isp (fm 2 ) 




0.617 


0.612 


0.610 


0.614 


< r 2 > s l (fnr 


') 


0.384 


0.547 


0.625 


0.563 


10 5 (C? 2 + 2C[ 3 ) [<§] 




-2.6 


-0.56 


0.55 


-0.71 


10 5 (C? 2 + 2C[ 3 ) [75] 




-3.3 


-0.55 


0.48 


-0.75 


10 5 (C[ 2 + 2C[ 3 ) [ 7 £] 




-0.55 


0.11 


0.33 


0.15 


10 5 C[ 3 [ 7 S S ] 




-0.56 


-0.02 


0.15 


0.03 


10 5 (C[ 2 + 4C[ 3 ) [<$] 




-4.1 


-0.27 


0.99 


-0.08 


10 5 C[ 2 [c5, 7 y 




-1.5 


-0.52 


0.26 


-0.78 


io 5 C[ 2 [c5, 7 f g ] 




-1.1 


-0.84 


0.12 


-1.3 


F (MeV) 




87.7 


63.5 


70.4 


71.0 


P (4) 




0.136 


0.230 


0.253 


0.254 






-0.083 


0.226 


0.059 


0.048 


F (4) 




0.169 


0.157 


0.153 


0.159 


P (6) 




0.051 


0.063 


0.067 


0.061 


^ 




0.736 


1.005 


1.129 


0.936 


/4 4) 




0.006 


-0.090 


-0.138 


-0.043 


/4 6) 




0.258 


0.085 


0.009 


0.107 


$ 




0.687 


0.938 


1.055 


0.874 


(4) 




0.007 


-0.100 


-0.149 


-0.057 


(6) 




0.306 


0.162 


0.094 


0.183 


«(2) 




0.734 


1.001 


1.124 


0.933 




-0.052 


-0.151 


-0.197 


-0.104 




0.318 


0.150 


0.073 


0.171 
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Figure 24: The result for the scalar radius < r 2 >^ as a function of L\ and Lg, with the 
dispersive estimate (abm) and the ChPT calculation. Figure from Ref. [142J. 



ml±/{ml±) QCD 

m W( m K±)QCD 



(2) (4) (6) 



(166) 



Notice that in |lU7j the corresponding numbers were quoted for the "Main Fit" which used the old K e4: 
data. In Ref. p.42] several of the order p 6 constants were also determined from the curvature in the 
form-factors, these are the combinations of C\ given in Tab. H3 The symbol in square brackets indicates 
the curvature of the form-factors used in the determination. The curvature is for the form-factor at 
zero normalized to 1, for the quantities labeled c. For the quantities labeled 7, the form- factors at zero 
are normalized to B or to its lowest order value. It is obvious from the table that there are some 
discrepancies still to be understood since not all determinations are compatible. The table shows the 
very large corrections to some of the scalar form-factors at zero. 

A glance at Table 01 shows that the pion decay constant in the chiral limit F can be substantially 
different from the value of about 87 MeV for the case of fit 10. The reason is that L\ is not that well 
constrained and thus allows for these different values. As discussed in the two following subsections, 
the constraints from im scattering and itK scattering could resolve this partly but the constraints are 
only marginally consistent. The consistent area does have values of F very close to the one of fit 10. 
A full conclusion cannot be drawn as a comprehensive analysis of all order p 6 constants is still lacking. 



4.9 Pion-pion scattering 

We have already treated pion-pion scattering in two-flavour ChPT. In order to complete the check of 
going from Ki& form-factors to pion-pion scattering, the latter process also needs to be known to order 
p e in three-flavour ChPT. This was accomplished in Ref. |151j . All pion-pion scattering in the isospin 



50 



limit can be described by the function A(s,t,u) defined in Eq. ([103)1 . 

This function was calculated in two-flavour ChPT to order p 6 in Refs. J72 EI] and is also known to 
order p 6 in three-flavour ChPT. The expression can be found in Ref. |151j . 

The convergence of the expansion is similar to the convergence in two-flavour ChPT. The various 
isospin amplitudes can be rewritten in terms of A(s, t, u) via 

T°(s,t) = 3A(s,t,u) + A(t,u,s) + A(u,s,t), 

T\s,t) = A(t,u,s)-A(u,s,t), 

T 2 {s,t) = A(t,u,s)+A(u,s,t), (167) 

where the kinematical variables t, u can be expressed in terms of s and cos 9 as 

£=--(s-4m*)(l-cos0), u = --(s-4m£)(l + cos0). (168) 

The amplitudes are expanded in partial waves using 

oo 

T\s, t) = 32vr J2(2Z + l)P/(cos 0)t\(s) . (169) 

1=0 

Near threshold these can be expanded in terms of the threshold parameters 

t\ = q 2i (a\ + bW + Otf)) , q 2 = \(s- Ami) ■ (170) 

Below the inelastic threshold the partial waves satisfy 



lmtl(s) = a(s) 



t\{s) 



\ a(s) = Jl-^. (171) 



In this regime, all partial waves can be written in terms of the phase-shifts. 

t\{s) = . l I {e 2 ^ - 1} . (172) 

y/l - (4ml/ s) 2l 

Some results of the phase-shifts were shown already for two-flavour QCD. 

The three-flavour calculation results for the threshold parameters are shown in Tab. Ejfor the sets 
of input parameters fit 10, A, B and C described in earlier subsections and in Tab. El 

It can already be seen from Tab. Hjthat Oq is very well predicted for sets of input parameters but 
Oq is somewhat smaller than the result of Ref. (HH]. Ref. |151j performed a full analysis, using the fits 
of |142j with the whole range of L r A , U & as input to determine the best fits. It was found that a° always 



fitted nicely but that the value of a^ restricted the values to a corner of the U A , L^ only. In Fig. 1231 the 
dispersive results from Ref. jHSl are shown together with the three-flavour order p 6 ChPT results. 

Similar results can be obtained for the other threshold parameters, but especially the higher threshold 
parameters have less convergence and have not been used as constraints. In Ref. |8*6 ^ll51j one also defined 
the subtreshold parameters C\ and C^. They can similarly be used to obtain constraints on the ChPT 
input parameters. The total constraints from pion-pion scattering on the input parameters are shown 
in Fig. |2i 
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Table 4: The values of the threshold parameters defined in ()170|) for the values of the input 
parameters of fits 10 |107j and A,B,C |142| . The lowest order values and the contributions 
from the order p 6 LECs, C[, are included. The threshold parameters are given in the 
corresponding power of m 2 + ■ Note that a\ and b\ are always given with the same power of 
ten. For fit 10, the three orders are quoted separately so the convergence can be judged. 
Table adapted from |151j . 







fit 10 


fit A 


fit B 


fit C 




p 2 


p4 


p 6 


total 


total 


total 


total 


0° 


0.159 


0.044 


0.016 


0.219 


0.220 


0.220 


0.221 


b°o 


0.182 


0.073 


0.025 


0.279 


0.282 


0.282 


0.282 


10 al 


-0.454 


0.030 


0.013 


-0.410 


-0.427 


-0.433 


-0.428 


10 fc 2 


-0.908 


0.151 


0.025 


-0.731 


-0.755 


-0.761 


-0.760 


10 a\ 


0.303 


0.052 


0.031 


0.385 


0.388 


0.389 


0.389 


10 b\ 


— 


0.029 


0.038 


0.067 


0.064 


0.063 


0.063 


\tfa\ 


— 


0.153 


0.080 


0.233 


0.223 


0.220 


0.221 


10 2 6° 


— 


-0.040 


0.007 


-0.033 


-0.035 


-0.036 


-0.036 


10 3 a 2 
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0.327 


-0.106 


0.221 
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0.221 


10 3 6 2 
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-0.385 
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io 4 4 
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0.20 


0.44 


0.64 


0.62 
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Figure 25: The scattering lengths Oq and a 2 , as a function of the input values L\ and Lg 
with the other L\ simultaneously refitted to the K e ± form-factors, (a) a° calculated to order 
p A and p 6 . The central value of jHH] of 0.220 is also shown, (b) a 2 , at order p 4 and p 6 . The 
two horizontal planes indicate the allowed region obtained in 86J. Figure from Ref. |151j . 
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Table 5: The results for the scattering lengths and ranges and the amplitude at the Cheng- 
Dashen point as well as the dispersive result. The scattering lengths and ranges are given 
in units of m n +. Table adapted from Ref. |157| . 





Fit 10 


HS2I 


1/2 

a 


0.220 


0.224 ±0.022 


10a 1 / 2 


0.18 


0.19 ±0.01 


in 3 / 2 
10a 


-0.47 


-0.448 ± 0.077 


10 2 a? /2 


0.31 


0.065 ± 0.044 


106 /2 


1.3 


0.85 ±0.04 


106 3 / 2 


-0.27 


-0.37 ±0.03 


1 CD 


2.11 


3.90 ±1.50 



4.10 Pion-kaon scattering 

An analysis of pion-kaon scattering similar to the combination of two-flavour ChPT at order p 6 and the 
Roy equations as done in Ref. [S3 EH] can in principle also be done for pion-kaon scattering. At the 
moment, this is not quite completed. The Roy equations analysis, extended to the pion-kaon scattering 
case using the so-called Roy-Steiner equations has been performed by the authors of Ref. |152j . Some 
earlier sum rule work and comparisons with ChPT can be found in Ref. |153j . 

On the ChPT side, the lowest order calculation was performed using current algebra methods in 
Ref. |154j . The order p 4 calculation can be found in Refs. |155[ 1156] and the full order p 6 calculation 
was performed by the authors of Ref. |157j . 

An alternative method using ChPT, is to treat the kaon as very heavy and only the pion as a 
Goldstone boson. This so-called heavy- kaon approach to pion-kaon scattering can be found in Ref. 158J. 

The pion-kaon scattering system can be decomposed into different isospin amplitudes and these in 
turn can be expanded in partial waves. The partial wave in turn are expanded around threshold and a 
full set of threshold parameters can be defined analogous to the pion-pion scattering case. There also 
exist subtreshold expansions in the pion-kaon case as well. These quantities have been calculated using 
the Roy-Steiner equations in Ref. |152j and to order p 6 in three-flavour ChPT in Ref. [157J. 

The entire comparison of inputs to ChPT, the resulting threshold and subtreshold parameters and 
their comparison with the dispersive results can be done as in the previous subsection. The results for 
the threshold parameters and thee inputs of fit 10 are shown in Tab. |S| As can be seen, good agreement 
exists for many of the threshold parameters. A full analysis is described in Ref. |157j and the constraints 
on the L4, Lg plane following are shown in Fig. |2B A region that is compatible with all results is 

L\ « 0.2 • 10~~ 3 and L\ « 0.0 • 10~ 3 . (173) 

A more extensive discussion can be found in Ref. |157j . 

One point deserving mention, the order p 6 calculation naturally obeys the chiral relations derived in 
the heavy kaon approach [158j . In particular, how to obtain the isospin odd low-energy theorem from 
the full order p 6 expressions of |151j is discussed in Ref. |159j . There has also been some recent progress 
in calculating the threshold parameters fully analytically |160j . 

4.11 7r,K->£wy 

To determine the parameter L\ at the same order as precision as was described in the previous sub- 
sections the processes n, K — > lv^ need to be calculated to order p e . This calculation has been done by 
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Figure 26: (a) The constraints on L\ and U % following from the dispersive results in tttt 
scattering jHHj as derived in |151j . (b) The same for the results from irK scattering [152 
as described in |157j . The curves from the boundaries are plotted on both plots to make 
comparisons easier. The shaded regions are excluded. Figure from Ref. |157| . 



the authors of Ref. |161j . They found a reasonably converging result. I refer to their paper for a longer 
discussion, but from Fig. E3 it can be seen that they found a nicely convergent series. 



4.12 Estimates of order p 6 constants 

One of the major problems is the sheer number of order p e constants that contributes. In some cases, 
they could be determined from experiment directly. Good cases here are those responsible for kine- 
matical quantities at order p 6 . These are the constants responsible for the curvature in the various 
form-factors, the kinematic factors 65 and b§ in pion-pion scattering and similar quantities in other 
decays. 

Especially the curvature in the vector form-factors leads to a well-determined combination. However, 
there are many more order p 6 LECs that cannot be so easily determined directly from data. Since at 
order p 6 typically many L\ contribute to particular processes, the determination of them becomes 
entangled between different processes. As a result, an estimate of an order p 6 LEC used in one process 
where an LI is determined sneaks in the determination of the other L\ and possibly order p 6 constants 
in the other processes. 

The solution, a full comprehensive analysis of all processes at the same time is a major undertaking 
which has not been done. It is also not clear whether sufficient data exist to constrain all involved L\ 
and order p 6 LECs. 

The solution which has mainly been used up to the present is to estimate the values of the order 
p 6 LECs by simple resonance saturation. This approach was started already in Ref. |2]. There it was 
shown that the experimental values of the order p 4 LECs cannot be the low-energy limit of the linear 
sigma model. They also found that exchange of vector mesons lead to reasonable predictions for the 
largest LECs at order p 4 . This work was systematized in Refs. |lb'2| where it was shown that the values 
of all order p 4 constants can be understood on the basis of one-resonance exchange. Ref. [163] showed 
how these results, when taken together with short- distance QCD constraints, are independent of the 
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Figure 27: The axial vector form-factor in ir — > iwy (left) and K — > £wy at order p 6 in ChPT. 
The dashed and solid curves represent the full results at order p 4 and order p 6 . Note that 
the latter includes data fitting. Figure from Ref. |161j . 



actual chiral representation chosen for the resonance fields. They also found that the main contribution 
came from the exchange of vector mesons. 

The approach used in most order p 6 papers consists in taking the simplest Lagrangian containing 
vector, axial- vector, scalar and sometimes pseudoscalar and tensor resonances, fitting its parameters as 
much as possible from data, and then taking the masses large to determine the LECs of ChPT. This is 
what has been done in basically all of the papers reviewed so far, with the exception of the two-flavour 
pion-pion scattering work. 



The main underlying theory argument is the large number of colours, N c , limit 



where QCD 



reduces to a theory of stable resonances. The additional assumption made in the estimates is that the 
lowest resonance in each channel will be important. 

Apart from the simplest version mainly used in the papers on ChPT to order p 6 more comprehensive 
approaches exist. There is the minimal hadronic approximation approach |165j . the ladder resummation 
model approach |166j and the resonance chiral model approach |lb'7t I146| I168| 1169] . All of these 
approaches are different ways of taking a certain number of resonances into account in each channel. 
They differ somewhat in the way the various couplings are treated. One problem is that there is 
an ambiguity in how one deals with short distance QCD constraints. It is known that there are 
incompatibilities for a finite number of resonances 



Other approaches are the older quark model approach, see |!7U| and references therein, as well those 
based on the Nambu-Jona-Lasinio model |171j . see |172] 11731 H7'4|. ITTo] and references therein. All these 
approaches give similar results in vector dominated channels but in those cases where quark mass effects 
play a major role, no full comparison with experimentally determined higher order coefficients has been 
done. For predictions regarding order p 4 LECs, most of these get similar results. 

One problem that affects all the large N c based methods is that the ChPT LECs are subtraction 
scale dependent while the estimates are not. This can only be solved by going beyond the leading terms 
in 1/N C also in the resonance models. Some attempts where also more references can be found are in 
Refs. fT7EllT771IT7S| . 

This area of estimating higher order p 6 LECs is where most improvement is needed in the future of 
higher order ChPT calculations. 
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5 Other 

In this section I dicuss the order p 6 calculations which have been done but which did not fit in the two 
main parts discussed in the previous two sections. 

5.1 Finite temperature and volume calculations 

The field of finite temperature calculations was started by Gasser and Leutwyler in Ref |179j . It means 
making the time part periodic and rotating it into the imaginary time direction. Two-loop calculations 
I am aware of in this area are the density dependence of the pion propagator by Schenk [180, as well 
as the pion mass and decay constant by Toublan |181j . 

At finite volume, there exists also work. Finite volume in ChPT was introduced by Gasser and 
Leutwyler in Ref. 1182J . A lot of work has been done to higher orders in this area using Luscher's 
method to determine the leading finite volume correction |183j . Proper two-loop calculations at finite 
volume in ChPT have only become available recently. The mass and the decay constant in two-flavour 
ChPT can be found in Ref. |184j and the quark-antiquark condensate in three-flavour ChPT in Ref. 185J. 

There is another regime at finite volume, the small volume or e regime. Here the lowest mode of zero 
momentum dominates the higher order corrections. It therefore needs to taken into account exactly and 
forms in this way a all-order calculation. This is a large area of active research and I will not discuss it 
here. 

5.2 Partially quenched calculations 

A large effort goes into numerically evaluating the functional integral of QCD. This is area is known 
as lattice QCD. One of the things that appears naturally in the way those calculations are done is 
the different treatment of quark lines connected to external legs, so-called valence quark lines, and the 
closed quark loops, called sea-quark loops. 

The effect of the latter is numerically very difficult to compute and has therefore often been approx- 
imated. Neglecting the sea-quarks completely is known as the quenched approximation and treating 
them with different masses from the valence quarks is known as the partially quenched approximation. 

In neither of these cases is the resulting theory a bona-fide field theory and the arguments that derive 
ChPT from QCD, given clearly in Ref. 0U] do not all hold. But both quenched and partially-quenched 
are still well defined statistical models and many of the arguments when thought of in terms of Feynman 
diagrams still seem to hold. We can thus hope that a ChPT extended to the case of quenched and 
partially quenched QCD still makes sense. Many of the qualitative prediction predicted by quenched 
ChPT (QChPT) and partially quenched ChPT (PQChPT) seem to be present in actual lattice QCD 
calculations. Especially the quenched chiral logarithms seem to exist. 

The extension of ChPT to the quenched case was started by Morel |186j and properly extended 
by Sharpe [1871 I188J and Bernard and Golterman |189j . Infinities and Lagrangians at one-loop have 
been discussed in Ref. 119(1]. The extension to the partially quenched case was again done by Bernard, 
Golterman and Sharpe |1911 I192J . Articles which provide a good overview of the field at order p A are 
Refs. [T331IT33J . 

Work on extending the work to two-loops started some years ago but there were two difficulties to 
be solved. One is to find the checks coming from the knowledge of the divergence structure and the 
Lagrangian at order p 6 and the other is the sheer complexity of the calculations in terms of the size 
of the expressions. The solution to the former turned out to be easy. Already in Ref. |19()j it was 
noted that the quenched Lagrangian and infinity structure at order p 4 had a very strong resemblance 
to the one of n^-flavour ChPT at the same order. More arguments for this were given in Ref. 195J. 
When Bijnens, Danielsson and Lahde started looking at two-loop PQChPT, they realized that a very 
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simple recipe allowed to obtain the full order p 6 both for the Lagrangian and the divergence structure. 
To understand this, one goes back to the suggestion of Morel |186j how to deal with the quenched 
approximation systematically. This is the approach worked out properly by Ref. |189j . In order to 
systematically remove the closed quark loops, one adds for each quark a bosonic quark, spin 1/2 but 
bosonic. Due to the different statistics, in this way all closed loops are exactly canceled. Partially 
quenched can be treated by leaving some of the quarks without a bosonic partner. Closed loops from 
these quarks then are not canceled. The formalism used in ChPT can be fully extended to this case by 
using supertraces instead of traces and adding extra bosonic and fermionic "Goldstone bosons." It was 
then also shown that in the partially quenched case the equivalent of the singlet eta would be heavy 
and could be systematically removed from the theory |194j . With that in hand, one realizes that all 
manipulations done in n^-flavour ChPT to obtain Lagrangians and divergence structures go through 
with the changes mentioned above and np changed to the number of active sea-quarks. The PQChPT 
Lagrangian at order p 6 and the divergence structure can thus be derived immediately from the work of 
Refs. |51| l53j. One also sees immediately the origin of the extra order p 4 term of |19bJ . The number of 
terms in the Lagrangian this leads to is given in Tab. El 

The complexity of the calculations was solved by combining the power of FORM |197j with a large 
amount of hand-optimized simplification of the various terms. At present no direct comparison with 
data are made, so I only show a representative plot. The relative corrections to lowest order for the 
case with three sea-quark flavours and valence-quarks of equal quark- mass, m 4 , and sea-quarks of equal 
quark- mass, 1712, are shown in Fig. EHJ Some visible features are the sizable corrections and the fact 
that the corrections do not vanish when mi —>■ but 777.4 kept fixed. The latter is the manifestation 
of the quenched chiral logarithm. More details can be found in the actual papers [198| I199"} I200[ f2 1 j . 
The actual calculations that have been performed are the masses for three sea-quarks |198[ 1201] , the 
decay constants for three sea-quarks |199| I201J and the same quantities for two sea-quark flavours as 
well 1200. The expressions can be found on the website [TU] . 

This is clearly an area where further progress is desirable. 



6 Conclusions 

In this review I have briefly discussed the foundations of ChPT and some of its technical aspects. The 
longest parts has been concerned with the actual calculations which have been performed at order p 6 
and their confrontation with experimental results. 

Here the two-flavour part, Sect. El is in rather good shape. Most calculations have been done and 
in many cases the contributions from unknown order p 6 constants are expected to be rather small. I 
particularly want to emphasize again here the work on pion-pion scattering where major progress in 
the theoretical description was made. 

In the three-flavour sector, Sect. EJ many calculations have been done and fitted to experimental 
results. One main highlight here is the predictions made in the context of K^ and its relation to the 
determination of V us but many other low-energy processes involving kaons and etas are also known. I 
have only briefly discussed the main remaining block towards a full comprehensive application of ChPT 
at NNLO, estimating in a consistent fashion all needed order p 6 LECs. Work in this area is going on, 
so I hope there will be significant progress in the not-too-distant future. 

And last, I have mentioned the finite temperature, volume and partially quenched work to two- loop 
order in ChPT. The interface between lattice QCD and ChPT is a growing area of research and should 
benefit both communities. At present only a few full two-loop order calculations in this sector are 
available as discussed in Sect. El 
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Figure 28: The relative corrections to the charged pseudoscalar meson mass, A M , and decay 
constant, Ap, to NNLO for the case of three sea-quarks. The valence quarks have quark- 
mass rrii and the sea-quarks mass m^. The plot is as a function of the quark- masses in the 
combination Xi — 2B mi and %4 — 2£> m 4 . The quantity plotted represents the sum of 
the NLO and NNLO corrections normalized to lowest order. The difference between two 
successive contour lines in the plots is 0.10. The values chosen for the LECs correspond to 
"fit 10" of Ref. [M]. Figure from Ref. J27TT] . 
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A A few remarks on alternative notations 

I have tried to stick in this review to one choice of notation. Many alternative parameterizations have 
been around especially for the two-flavour case. 

There are two main conventions around for the parameters F, F , F which differ by a factor of y/2. 
Which is in use can normally be detected by the forefactor in Eq. (|48|) . F 2 /A corresponds to an F n of 
order 93 MeV, F 2 /8 to an F n of order 130 MeV while in the older literature sometimes even an F n of 
order 180 MeV appears. 

The parameters B, Bo, B often appear in the combination F 2 B/2 which is frequently called r. 

Many people use a left-right version instead of the right-left version introduced in Ref. [3]. They 
typically use S = IP as the quantity containing the Goldstone bosons. The quantity u also appears 
often denoted as £ and exists as well in a version with left and right interchanged. 

For the two-flavour case, the fact that 

SU{2) x SU(2)/Z 2 = SO{4) (174) 

and 

SU(2)/Z 2 = SO{3), (175) 

allows a parameterization inspired by the breaking of SO (4) to 5*0(3). These are sometimes called 
sigma-model parameterizations. In Ref. [2j a four-vector U = (U° , U 1 , U 2 , U 3 ) was used with the 
restriction U T U = 1. Here the pion fields appear via U l = it 1 / F for i = 1,2,3. All needed external 
fields can also be brought into this representation. The explicit connection between the exponential 
and the sigma-model representation is 

U l = ~(rV) (176) 

where the r % are the three Pauli-matrices. Weinberg often uses a three-vector fr, and a covariant 
derivative which depends nonlinearly on the fr field. This is the form he used in his original paper |37j . 
All of these parameterizations are equivalent, as it has been proven that all can be brought into one 
standard form in Ref. [SB]- However, when performing calculations, individual Feynman diagrams can 
be quite different in the different formulations. As an example, the pion wave function renormalization 
vanishes at NLO in the parameterization of Ref. |2] but it doesn't in the exponential parameterization 
used in this paper. 
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